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ABSTRACT 

H I line widths are typically interpreted as a measure of ISM turbulence, which is 
potentially driven by star formation. In an effort to better understand the possible con- 
nections between line widths and star formation, we have characterized H I kinematics in 
a sample of nearby dwarf galaxies by co-adding line-of-sight spectra after removing the 
rotational velocity to produce an average, global H i line profile. These "superprofiles" 
are composed of a central narrow peak (~ 6 — 10 km s~^) with higher- velocity wings 
to either side that contain ~ 10 — 15% of the total flux. The superprofiles are all very 
similar, indicating a universal global H I profile for dwarf galaxies. We compare charac- 
teristics of the superprofile to various galaxy properties, such as mass and measures of 
star formation (SF), with the assumption that the superprofile represents a turbulent 
peak with energetic wings to either side. We use these quantities to derive average scale 
heights for the sample galaxies. When comparing to physical properties, we find that 
the velocity dispersion of the central peak is correlated with (Shi). The fraction of 
mass and characteristic velocity of the high velocity wings are correlated with measures 
of SF, consistent with the picture that SF drives surrounding Hi to higher velocities. 
While gravitational instabilities provide too little energy, the SF in the sample galaxies 
does provide enough energy through supernova, with realistic estimates of the coupling 
efficiency, to produce the observed superprofiles. 
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Introduction 



The neutral hydrogen (Hi) component of the interstellar medium (ISM) is an ideal tracer 
of kinematics in disk galaxies. Hi is observable in some galaxy disks far beyond optical emission, 
making it a superb tool for probing large-scale kinematics speeds. On smaller scales, the Hi velocity 
dispersions offer a way to measure the random turbulence velocities on scales of 10 — 200 pc. By 
connecting the H I velocity dispersion with possible drivers of turbulence, we can study the sources 
of energy in the ISM. 

H I velocity dispersions typically vary between 5 — 15 km s~^ across a wide range of disk galaxy 
types, and generally decrease in the outskirts of galaxies to ~ 6 — 10 km s~^ (e.g., Tamburro et al. 



2009 ) . Generally, these line widths are thought to be due to turbulent velocities rather than thermal 
broadening. For example, Wolfire et al. (1995) found two stable temperatures for Hi gas: ~ 150 K 
for the cold neutral medium (CNM) and ~ 7000 K for the warm neutral medium (WNM). These 
temperatures correspond to velocity dispersions of ~ 1 km s~^ and ~ 7 km s^^ at typical ISM 
pressures, which is often smaller than the observed line widths in nearby galaxies. This mismatch 
suggests that the line widths are set primarily by turbulence. However, the time scale for dissipating 
turbulent energy is ~ 10^ yr (Mac Low 1999). Energy must therefore be continually injected in 
order to maintain the Hi line widths we see in galaxies. 

The sources of energy that drive turbulence are still debated. A number of studies have 
suggested that star formation can provide the necessary energy to generate Hi turbulence in the 
inner regions of galaxies (e.g., Kim et aL]|1998 Tamburro et al.||2009 Joung et al.||2009 ). However, 



Hi velocity dispersions are still substantial at large radii, whereas the majority of star formation in 
galaxies is contained within r25, defined as the radius where the B-band surface brightness drops 
below 25 mag arcsec"^ (e.g. 



Kennicutt 


1989 


Bigiel et al. 


2010 



rate (SFR) falls off much more quickly than the H l velocity dispersion, implying that it cannot be 
the only contribution to the Hi velocity dispersion in disk galaxies. 



Other proposed drivers of turbulence are the magneto-rotational instability (MRI; ISellwood 



&: Balbus 1999), shear from rotation curves (e.g., Schaye 2004), or gravitational instabilities (e.g.. 



Wada et al. 2002). One can potentially distinguish among these various mechanisms by compar- 
ing the observed energy in turbulence to the energy available from the possible drivers. Many of 
these processes should be effective in spiral galaxies, which have spiral arms and exhibit differen- 
tial rotation, but should be less strong in their lower mass dwarf counterparts, which lack spiral 
structure and show solid-body rotation. However, the observed Hi velocity dispersions of spirals 
are surprisingly similar to those of dwarfs. 

In this paper we take a different approach to study the global behavior of the Hi velocity 
dispersion and its relationship to possible drivers of turbulence. By working on global scales, we 
can not only increase the signal-to-noise of individual spatially-resolved line-of-sight spectra but 
also remove the assumption that input energy must necessarily couple to Hi in the same spatial 
region. Our work extends the many previous studies of turbulence in the ISM but uses better 
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data over a wider baseline in galaxy mass. We also improve the characterization of the average 
velocity dispersion. In contrast, most earlier papers typically use the intensity-weighted second 



(e.g., 


Tamburro et al. 


2009 


Petric & Rupen 


2007 


), or I 



(e.g., van Zee &: Bryant 1999). However, the second moment can be artificially increased by gas 
with anomalous velocities, such as bulk inward or outward flows or expanding Hi, while single 
Gaussian fits are unable to represent asymmetric line-of-sight spectra. In addition, most literature 
studies of H I turbulence have not used a uniform sample of observations, as few such samples have 



been available until recently. Instead, studies focused on a single galaxies (e.g., Petric Sz Rupen 



2007) or worked at the instrumental resolution for each galaxy (e.g., Tamburro et al. 2009). The 
combination of such studies means that H I turbulence is sampled on different physical scales in each 
galaxy. Since turbulence is larger on larger physical scales (e.g., Zhang et al.||2012 ), this mismatch 
in spatial resolution makes galaxy-to-galaxy comparisons dubious. 

A better, uniform measurement of the typical underlying Hi turbulence is necessary to ac- 
curately constrain the detailed kinematics of the ISM. By co-adding Hi line-of-sight spectra after 
removal of the rotational velocity, we can obtain an average measurement of H l turbulent velocities. 
A small number of previous studies have followed a similar approach as we undertake here. |Dickey| 



et al. (1990) found relatively constant Gaussian line widths in the face-on spiral NGC 1058, with 
median profiles at some radii exhibiting wings larger than expected from a simple Gaussian profile. 
These high- velocity wings were shown to exist in the average line profiles regardless of the average 
FWHM of the contributing line-of-sight profiles ( Petric Sz Rupen]|2007 ). Similar results were found 



by Boulanger fc Viallefond (1992) and Kamphuis k Sancisi (1993) in NGC 6946. Braun (1997) 



also found H I gas at higher velocities compared to the average central H i line width in a number 
of other spirals by studying average Hi line profile shapes. However, the existence of high- velocity 
wings superimposed on a Gaussian center may not be ubiquitous; no evidence of such wings is seen 



in NGC 5457 ( |Rownd et al.||l994D or the outer regions of NGC 1232 ( [van Zee Bryant||l999[ ). 
Unfortunately, the majority of the studies of average Hi line profiles had poor velocity resolution 
(> 5.2 km s-^ Braun 1997 ), coarse spatial resolution (> 1 kpc; [Dickey et al.||l990[ [Petric &: 



Rupen[[2007t [van Zee fc Bryant||1999[ ), or both ( [Boulanger Viallefond]|1992[ [Rownd et aL][1994[ ) 



A recent study by lanjamasimanana et al. (2012 hereafter 12012), generated average Hi profiles 
( "superprofiles" ) using a similar approach as this paper, but for a number of more massive spirals 
within D ~ 10 Mpc. They found the same basic line profile structure as we see here, and proposed 
that they may be comprised of emission from the CNM and WNM. 

Recently, a number of Hi synthesis observation surveys of nearby galaxies have greatly im- 
proved the available data. Compared with the numerous published Hi studies of single galaxies, 
these surveys can provide better spatial and spectral resolution as well as a uniform observing 
setup. The Hi Nearby Galaxy Survey (THINGS; Walter et al. 2008) pioneered this new era of 
high-resolution Hi surveys by observing 34 nearby spiral galaxies with high spatial (6 — 10") and 
spectral (1.3 — 5.2 km s~^) resolution. The Very Large Array ACS Nearby Galaxy Survey Treasury 
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Project (VL A- ANGST; |Ott et al.|[2012[ ) followed in its footsteps, extending THINGS to smaller 
galaxy masses at a similar sensitivity and better spectral resolution (0.6 — 2.6 km s~^). Other sur- 



veys, such as FIGGS ("Faint Irregular Galaxy GMRT Survey" ; [Begum et al.||2008[ ) and LITTLE 
THINGS ("Local Irregulars That Trace Luminosity Extremes-THINGS" ; [Hunter et al.|[20T2l ) have 
also sought to provide a uniform sample of H i observations of dwarfs with similar observing setups. 

In this paper we focus on the global measurements of H i kinematics in a wide range of dwarf 
galaxies chosen from VLA-ANGST and THINGS. We present a method to measure the intrinsic 
H I kinematics in low-mass disk galaxies by co-adding flux- weighted H i line profiles after removal of 
the rotational velocity. The combined sample covers a wide range of galaxy properties, allowing us 
to examine the overall H i kinematics in a broader range of environments than previously studied. 
In § [2| we describe the data we use to determine galaxy global properties. We next explain our 
method of characterizing the global H i gas kinematics in § [3[ We then discuss our parameterization 
and the physical interpretation of these superprofiles in § [4[ In § [5| we investigate significant 
correlations between superprofile parameters and galaxy physical properties. In § [6| we discuss 
the potential physical cause behind correlations with each parameter; give energy estimates for 
driving kinematics in the different components; assess limits on Hi scale heights; and examine the 
possibility of a universal H i velocity profile for dwarfs. Finally, we summarize our results in § [7| . 



2. Sample and Data 



In the following sections, we briefly describe the Hi data available from VLA-ANGST and 
THINGS as well as our sample selection criteria. We then discuss our conversion from the data 
to the physical properties of our sample galaxies. If available, we use the ANGST TRGB distance 



from Dalcanton et al. (2009). Otherwise, we use the distances compiled in Karachentsev et al 
(2004). When necessary, we correct all published quantities to our adopted distances, and we 
include published or estimated distance uncertainties in the uncertainties for all our calculated 
quantities. In § [2.1.3 we discuss the criteria that we use to select the galaxies for our sample. 



The 9 galaxies from THINGS and 14 galaxies from VLA-ANGST in our sample are listed in 
Table [T] in order of decreasing baryonic mass along with their basic physical properties. We list 
(1) the galaxy name; (2) alternate names; (3) Hi survey (VLA-ANGST or THINGS); (4-5) right 



ascension and declination; (6) distance; (7) inclination i, (see § 2.3.1); (8) total Hi mass taken from 
Walter et ah] ( |2008[ ) or |Ott et al.] ( |2012[ ); (9) r25, from [Karachentsev et ar] ( [2004[ ); (10) inclination- 



corrected width at 20% of the total line profile, W20, from Walter et al. (2008) or Ott et al. (2012) 



and corrected using i as listed in this table; and (11) de Vaucouleurs T-type taken from Walter 



et al. ( 


2008) or 


Ott et al. 


(2012 



If necessary, previously-published quantities such as mass are 



corrected for the adopted distance. 
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2.1. HI Data 

We use a combination of NRAO Very Large Array (VLA) H i data from both THINGS and 
VLA-ANGST for our analysis. The two surveys provide complementary information, as THINGS 
targets are primarily large spiral galaxies while VLA-ANGST probes gas-rich galaxies at lower 
mass scales. Our final sample is composed primarily of dwarfs that span a wide range of galaxy 
properties, including absolute magnitude; stellar, gas, and baryonic masses; SFRs; and rotation 
speeds. 



We use the robust-weighted data cubes published in Ott et al. (2012) and Walter et al. (2008) 
from VLA-ANGST and THINGS, respectively. This weighting scheme offers ~30% better spatial 
resolution and a well-behaved synthesized beam with only a moderate decrease in sensitivity. Both 
VLA-ANGST and THINGS provide two sets of robust-weighted data cubes, the standard and flux- 
rescaled cubes. Standard data cubes have uniform noise properties but incorrect fluxes, while the 
flux-rescaled cube has been scaled to have correct fluxes in exchange for more complicated noise 
properties. Because we use both cubes in our analysis, we specify which cube we are using in each 
step. To ensure that we are sampling the same physical scales in the galaxies' ISM, we work at a 
common physical resolution of 200 pc as discussed in §2.2[ 



The parameters of the Hi observations are listed in Table [2j We list (1) the galaxy name; (2) 
velocity resolution, At;; (3) FWHM beam corresponding to 200 pc physical resolution, 02OOpc; and 
(4) rms noise in a single channel of the 200 pc convolved standard cube, cTchan- 



2.1.1. VLA-ANGST 



The ACS Nearby Galaxy Survey Treasury (ANGST) Program obtained multi-color HST pho- 
tometry of a volume limited sample of galaxies within 4 Mpc, excluding the Local Group, and 



provides spatially-resolved star formation histories for its sample (Dalcanton et al. 2009). As a 



followup, the VLA-ANGST survey targeted all galaxies in ANGST that were visible with the VLA 
{6 > —30°), showed signs of having observable Hi reservoirs, and lacked adequate previous Hi 
observations. Many of the galaxies in the sample of 35 are therefore low-mass, low-luminosity 
dwarfs. 

The velocity resolution of the survey is 0.6 — 2.6 km s~^, which is necessary to study the 
detailed H i line profiles in galaxies with the low peak rotation speeds characteristic of the sample 



(e.g.. Warren et al. 2012 ). The typical instrumental spatial resolutions at the median galaxy distance 
is ~ 7", corresponding to ~ 100 pc at the median distance of 2.8 Mpc. 
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2.1.2. THINGS 

THINGS provides a complementary sample of 34 large, gas-rich spirals chosen mainly from the 
Spitzer Infrared Nearby Galaxy Survey (SINGS). Since THINGS galaxies have characteristically 
higher masses and rotation speeds, coarser velocity resolutions were often required to fully cover 
the Hi emission. Therefore, the velocity resolution is often coarser than VLA- ANGST, and the 
majority of THINGS observations have velocity resolutions of either 2.6 or 5.2 km s~^. THINGS 
resolves spatial scales of ~ 7", corresponding to ~ 200 pc at the median distance of ~ 6 Mpc. 



2.1.3. Sample Selection 

To increase the robustness of our results, we select a high-quality subset of the 63 detected 
galaxies in VLA- ANGST and THINGS surveys for our analysis. We consider only disk-dominated 
galaxies (de Vaucouleurs T-type > 3) to avoid confusion with the bulge; all selected galaxies are 
7 < T < 10. 

We further exclude galaxies that suffer from one or more of the following problems, with the 
number of galaxies eliminated due to each criterion given in parentheses. 



1. Instrumental physical resolution larger than our working resolution of 200 pc (18 galaxies). 

2. Velocity resolution Av > 5.2 km s~^, which complicates determination of the peak velocity 
and approaches the width of turbulent regions in these galaxies (16 galaxies). 

3. Inclinations > 70°, which could lead to artificially broadened line profiles due to beam smear- 
ing (7 galaxies). 

4. Noticeable contamination from the Milky Way or from a companion, which would hinder 
separation of the galaxy Hi emission from its companion (8 galaxies). 

5. Fewer than 10 independent beams above the signal-to-noise threshold where we can accurately 
measure Upeak (S/N > 5; see § 3.2 ) at our working resolution. Galaxies with fewer independent 
beams show very noisy co-added profiles and have more than 50% of H i flux in pixels below 
our S/N threshold (8 galaxies). 

6. A lack of ancillary far-ultraviolet (FUV) imaging, needed to determine the current average 
star formation rate (10 galaxies). 



These cuts eliminate 40 potential galaxies from our sample; the majority of these were cut 
because they failed item (1) or (2) of the above criteria, and a number of galaxies failed more than 
one of the criteria. The final sample for analysis has 23 galaxies, with 14 from VLA-ANGST and 
9 from THINGS. 
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2.2. Hi Data Preparation 

To provide the best galaxy-to-galaxy comparison of H i kinematics, we must take the spatial 
resolution into account. Since the velocity dispersion is typically larger on larger spatial scales (e.g.. 



Zhang et al. 2012), we must ensure that we are sampling the Hi kinematics on the same spatial 
scale in our sample galaxies. Since the instrumental angular resolution is roughly the same but the 
galaxies are at different distances, we must apply spatial smoothing to some of the cubes to ensure 
that the same spatial resolution is sampled from galaxy to galaxy. This is essentially equivalent to 
placing all galaxies at the same distance. 



We choose a spatial resolution of 200 pc to match that used by Warren et al. (2012), allowing 
for future comparison of our results. This resolution is also a good compromise between potential 
sample size and physical resolution at the distance of each galaxy, which is limited primarily by 



galaxy distance. It also matches results from Joung et al. (2009), who find that most turbulence is 
contained on spatial scales of 200 pc or less. 

To apply the spatial smoothing to the data cubes, we first calculate the beam size that yields 
200 pc resolution at each galaxy's distance. We produce spatially-smoothed versions of both the 
standard and flux-rescaled data cubes at 200 pc resolution using the AIPS task CONVL, which 
accounts for the original beam major axis, minor axis, and position angle. 

To generate a mask for these 200 pc data cubes, we first convolve the original standard cube 
to 45" resolution (convl). Next, we measure a , .^n, the rms noise in the 45" cube, and mask 

^ ' ' cnan,45 ' ' 

all emission below 2>a , .ui using the AIPS task blank. Finally, we remove any remaining non- 
emission regions by hand. To regenerate the moment maps for the 200 pc resolution data, we 
blank the convolved, flux-rescaled data cube in regions outside of the mask. We then use the AIPS 
task XMOM to produce zeroth, first, and second moment maps. For the remainder of the paper, 
all mention of data cubes or moment maps refer to the 200 pc data sets described in this section, 
unless otherwise specified. 

We note that convolution to a circular beam means that inclined galaxies have slightly larger 
physical resolution along their minor axis than along their major axis. However, the uncertainty 
of inclinations and position angles for the majority of galaxies in our sample makes it difficult to 
correct. Therefore, we choose to use the simplest option of a circular beam. 

Finally, we include in our final analysis only line-of-sight spectra with a signal-to- noise (S/N) 
> 5, where S/N is defined as the ratio between our fits to the peak divided by the rms noise in the 



line-free channels of the data cube. The reasons for this choice are discussed further in § 3.2 



- 8 - 



2.3. Converting Data to Physical Quantities 



In this section we describe the methods we use to measure the physical quantities discussed 
in the paper so that we are able to compare them to the H i superprofile properties of the sample. 
For many of the quantities, such as the star formation rate (SFR) or the Hi surface density ShI; 
we calculate the global properties using only the pixels whose Hi line-of-sight spectra contribute 
to the superprofile (S/N> 5). This choice provides a matched aperture measurement that allows 
us to consider only regions that are able to directly affect the H i measured by the superprofiles. A 
notable exception is the total baryonic mass of the galaxy, Mbaryon,tot) which we use as a proxy for 



halo mass (§ 2.3.2). In this case, including only pixels above the S/N threshold would artificially 
underestimate the total baryonic mass, and therefore the halo mass, of galaxies in the low signal- 
to-noise regime. The halo mass is not expected to directly influence Hi velocity dispersions, but it 
is useful to first characterize how the H i superprofile properties behave as a function of total halo 
mass before exploring their connection with other physical properties. 

We list our derived quantities in Table [Sj We give these quantities, followed by the relevant 

(3) Hi mass, Mhi, 



section: (1) the galaxy name; (2) Mb, 



aryorijtot ; 



M^, § ^J:7\ (5) SFR, § [2X5| (6) SFR / Mm 



2.3.2 



2.3.3, (4) stellar mass, 



2.3.6, (7) average star formation rate surface 
density, (Ssfr), § 2.3.5, and (8) (Shi), § 2.3.3 We measure the quantities listed in columns 3 — 7 



using only pixels above our S/N threshold. 



2.3.1. Galaxy Inclination 



For disk galaxies, the best inclination is usually the one inferred from tilted ring model fits 
to the Hi velocity field. Inclinations derived this way are available for 12 galaxies in our sample 



from a variety of sources in the literature (e.g.. 


Skillman et al. 


1!)SS 


Begum et al. 


2005 


de Blok 


et al. 


2008 


S waters et al. 


2009 


Oh et al. 2011); we use these inclinations if available. We note 



that for strictly solid body rotation, as is common in dwarfs, the velocity field does not contain any 
information about the inclination angle of the disk. 



However, 11 galaxies in our sample do not have previously-derived rotation curves and are 
therefore lacking these inclination estimates. In the absence of velocity field analysis, the traditional 
method is to measure ellipticity from B-band observations and then to derive an inclination after 
assuming an intrinsic disk thickness. This method often fails dramatically in dwarf galaxies, since 
it relies on the assumption that the intrinsic disk structure is well-traced by the B-band surface 
brightness. Unfortunately, SFRs in dwarf galaxies are much lower than in spiral galaxies, such that 
the massive stars that dominate the B-band surface brightness are formed stochastically across the 
disk. Since their light does not smoothly trace the galactic disk, we must turn to another indicator 
to measure the projected galactic disk. 

As opposed to B-band observations, near-infrared observations are dominated by flux from the 
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older stars and should therefore provide a better measurement of the projected shape of a galaxy's 
disk. All sample galaxies without inclinations derived from rotation curves are part of the Local 



Volume Legacy survey (LVL; Dale et al. 2009), which provides photometric infrared observations 
of galaxies within 11 Mpc. For these galaxies, we fit ellipses to isophotes in the the LVL 3.6/im 
Spitzer images. We then calculate inclination by assuming: 



sin^ i 



1 - jb/a)' 



(1) 



where (6/a) is the measured axial ratio and qo is the intrinsic disk thickness. We use the values for go 
provided in Karachentsev et al. (2004) for different galaxy types, which range between 0.12 — 0.2. 



The derived inclinations only change by < 5° when we increase qq to a fixed value of 0.3. We 
estimate an uncertainty Ab/a ~ 0.05 from repeated measurements of "best-fit by eye" ellipses to 
the 3.6/im surface brightness distribution, which leads to uncertainties in the inclination of ~ 5°. 
When we compare inclinations derived using this method to those of galaxies with tilted ring 
inclinations, we find that the inclinations typically differ by less than 10°. The inclinations derived 
using Hi morphology ( Begum et al.|[2008 ) for the five galaxies that overlap both samples are within 
5°, with the exception of DDO 187. We therefore estimate our total uncertainty on the inclination 



as (Tj 



10° 



We denote the galaxies whose inclinations have been derived using this method with * in 
column 7 of Table [TJ 



2.3.2. Halo Mass 



In large spiral galaxies, the inclination-corrected widths of the Hi integrated line profile at 
20% or 50% of the peak {w2o and w^o) are good tracers of the total halo mass, as the rotation curve 



flattens to approximately the circular velocity (e.g., Verheijen 2001). In low-mass dwarf galaxies, 
however, rotation curves often continue rising past the extent of the observable H i, so any measured 
velocity provides only a lower limit on the circular velocity of the halo. Second, the global profiles 
are generally Gaussian and, due to the small rotational velocities, are more affected by turbulent 
motions (Begum et al. 2006). The ability to derive the intrinsic velocity width also requires an 



accurate knowledge of the galaxy inclination, which is uncertain for many of the low-mass dwarfs 
in our sample. 



However, detailed studies of the baryonic Tully-Fisher relation in dwarf galaxies (e.g., Geha 
et al.|2006 Stark et al.|2009 ) indicate a strong correlation between baryonic mass, as measured with 
gas and stars, and halo mass, as measured with either ■W20 or w^q when inclinations are well-known. 
This correlation is expected if all halos in our sample have approximately the same baryon fraction. 
In the absence of reliable inclination-corrected W20 measurements, however, it is preferable to use 
the baryonic mass as a proxy for halo mass in dwarf galaxies. We therefore use -Mi5a,ryoii,tot 

to 

indirectly measure halo mass in our sample. 
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Unlike the other measurements described in this section, we use the total baryonic mass instead 
of only the mass contained in the same pixels as we use to derive the superprofiles. Because this 
measurement is simply a proxy for halo mass and does not have a direct causal connection to the 
Hi properties, we include the entire baryonic mass for each galaxy. Including only pixels above our 
S/N threshold would underestimate the baryonic mass, and therefore the halo mass, of galaxies in 
the low signal-to-noise regime. 



We calculate Mb; 



.aryon.tot = l-36MHi,tot + M^.tot, using listed values for MHi,tot from [Walter 
et al.l (120081) and lOtt et al.l (|2012|) ; and L3.6 values from ID ale et al.l (|2009|) and then apply the 



L3.6 — )• conversion in §2.3.7 The factor of 1.36 accounts for helium; we neglect metals as the 
low metallicities of our sample imply that < 1% of the gas is composed of metals (see § 2.3.4). In 
most cases, -/Vfbaryon,tot IS higher than the aperture-matched measurement of Mbaryom primarily due 
to Hi in low S/N pixels that do not contribute to the superprofiles. 



2.3.3. Hi Mass 

We calculate the average Hi surface density, (Shi), and Hi mass, Mhi, from the convolved Hi 
total intensity maps. We first convert this map to a de-projected surface density in pc~^: 

E„ (Me PC-) = 12.14 ,^H,i_^%HM„.„ 

where S^v is the Hi surface brightness in Jy beam~^ km s~^, FWHMmaj and FWHMmin are the 
beam major and minor axes in arcsec, and i is the inclination. 

To calculate the average (Sri) for each galaxy, we average the Shi map using all the pixels 



above our S/N threshold (see § 3.2). We then calculate Mhi by summing Shi times the de-projected 
physical area of each pixel. We note that the mass we calculate in this step is less than the total 
Hi mass of each galaxy given in Table [T] due to our S/N threshold. 



2.3.4. Gas Mass 

We calculate the gas mass by assuming Mgas = 1.36 Mhi. We include the factor of 1.36 to 
account for the presence of helium. As most of our galaxies are dwarfs with currently-undetectable 
H2, we neglect possible contributions from molecular gas. Only two galaxies in our sample have 
detected H2, and in both cases the molecular gas contributes < 10% to the gas mass. We also 
do not correct for metals, since the metallicities of our sample are likely to be low; the 4.5/im 



luminosity-metallicity relations given in Berg et al. (2012) imply metallicities of 12 + log (0/H) < 8 



for all but two galaxies (NGC 4214 and NGC 7793). Corrections for the masses of heavy elements in 
the ISM are therefore less than 1% for all of our galaxies, and thus much less than the uncertainties 
in the gas masses. 
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2.3.5. FUV + 24nm Star Formation Rate 



We measure the SFR by combining GALEX FUV and Spitzer 24/im luminosities following the 



prescription from Leroy et al. (2008). FUV emission primarily traces unobscured star formation 
that has occurred within the past ~ 10 — 100 Myr, while the 24/im emission traces warm gas that 
has been heated by embedded star formation on timescales of 3 — 10 Myr (Calzetti et al. 2007). 



Therefore, the combination of FUV with 24/um provides a measurement of both the embedded and 
the unobscured recent star formation. 



We generate pixel-by-pixel maps of SFR following the formalism described in Leroy et al. 



(2008). The empirically-calibrated relationship between FUV emission, 24^m emission, and SFR 



is given by: 



where SgpR is in Mq kpc 
conversion assumes a 



SsFR 
2 



(8.1 X 10^^/fuv + 3.2 X 10^^/24) cosi, 



(3) 



yr and the FUV and 24/xm intensities are in MJy ster . The 



Kroupa (2001) IMF with a maximum mass of 120 Mq as implemented in 



STARBURST99 (Leitherer et al. 1999). The measured SFRs are smaller by a factor of 1.59 com- 
pared to a Salpeter (1955) IMF with a mass range of 0.1 — 100 M0 when normalized for the same 
number of ionizing photons. At EsFR < 10~^ Mq yr"^ kpc~^ the 24;um emission is an upper limit 
to the SFR due to the diffuse dust component of the disk. However, only two of our galaxies have 
average (Esfr) < 10~^ Mq yr~^ kpc~^, and in these cases, over 90% of the SFR is determined by 
the FUV component alone. The 24//m tracer contributes less than 20% of the total SFR in all but 
four galaxies (DDO 53, NGC 2366, NGC 4214, and NGC 7793). 



We start with publicly available FUV and 24|Um images from LVL (Dale et al. 2009). These 
maps have resolutions of ~ 5" and 1.6" , respectively. For the FUV images, we subtract a small 
sky background and correct for Galactic extinction using the dust maps from Schlegel et al. (1998) 
and Ap\jy/E(B — V) = 8.376. We mask foreground stars, identified by pixels that have NUV / 
FUV flux ratios > 15. The 24/im maps have already been background-subtracted, so no additional 
correction is applied. We then convolve both the 24^m and FUV data to our working resolution, 
place both maps on the astrometric grid defined by the H i data, and apply Equation [3j 

To calculate the global SFR for each galaxy, we sum the Ssfr contribution from only the 
pixels above our Hi S/N threshold multiplied by the de-projected physical area of each pixel. To 
estimate (Ssfr), we divide the SFR by the de-projected physical area covered by pixels above our 
S/N threshold. 

We have attempted to calculate Ha-based star formation rates using LVL Ha maps, as they 
trace star formation on timescales ~ 10 Myr, shorter than the FUV+24/im measurement. However, 
the SsFR values implied by Ha observations are often below the 10~^ M0 yr~^ kpc~^ level where 



Ha maps are no longer reliable (e.g., Leroy et al. 2012). Additionally, at the low SFR typical of 



our sample, FUV tracers of SFR appear to be more robust than Ha tracers (e.g., Lee et al. 2011 



Leroy et al. 2012). We therefore use the FUV+24//m as our only SFR tracer. 
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2.3.6. Star Formation Rate per M}ii 



The star formation rate per unit Mgas is often taken to be the star formation efficiency (SFE). 
However, such an interpretation can be problematic in dwarf galaxies. The H2 component of the to- 
tal gas mass is notoriously difficult to measure in dwarfs but must exist if our current understanding 



of star formation is correct (although see Krumholz 2012). Using dust as a proxy for H2, Bolatto 



et al. (2011) find molecular gas fractions of Ehj/Shi 0.1 on 200 pc scales in the SMC. However, 



the typical Kennicutt-Schmidt relation for standard spirals tends to overpredict the star formation 



rate for a given total gas mass even on a pixel-by-pixel basis (e.g., Leroy et al. 2008 Bolatto et al 



2011). Such observations could be explained if dwarf galaxies have a fundamentally different SFE, 



or if the Hi surface density in dwarfs is simply less directly connected to star formation than in 
larger spirals. Since the Ssfr-5^H2 relation is comparable to that found in larger spirals, the latter 
is likely the case. Therefore, our SFR / Mhi measurement more likely traces the ability of recent 
star formation to affect the Hi gas than any true SFE effects. 



To calculate the global average (SFR/ Mhi), we simply calculate: 

SFR 

(SFR/Mhi) 



Mhi 



(4) 



where SFR and Mhi are the star formation rate and H i mass from contributing pixels as derived 
in 



2.3.5 and S 2.3.3 



2.3.7. Stellar Mass 



We use 3.6^m Spitzer data from LVL to estimate the stellar mass, using the method in Ap- 



pendix C of Leroy et al. (2008 ). This band primarily traces the light from older stellar populations. 



In more massive galaxies, the 3.6/im band can also contain emission from hot dust and polycyclic 
aromatic hydrocarbons (PAHs). However, dwarf galaxies show reduced PAH emission compared to 
larger galaxies of the same color (e.g., |Hogg et al.|2005t|Engelbracht et al.|2005t [Madden et al.|2006| 



Rosenberg et al. 2006; Jackson et al.|[2006 ). Since our sample is primarily composed of dwarfs, we 
do not account for PAH emission when converting from 3.6/xm intensity to stellar mass. 



We use the empirically-derived conversion from Leroy et al. (2008): 

Ik 



= T^,k(t^)/3.6 cos? 
-'3.6 



280 13.6 cosz. 



(5) 

/3.6 is the 3.6/im 



where T^^k ~ 0.5 is the mass-to- light ratio in the K-band, is in Mq pc ^ 
intensity in MJy ster~^ and {Ik/Is-g) ~ 1-81 is the 3.6//m-to-K-band conversion as derived by Leroy 



et al. (2008). The conversion assumes a Kroupa (2001) IMF. The mass-to-light ratio has a scatter 



of 0.1 dex. Further discussion of the conversion is given in Leroy et al. (2008). 



We start with the LVL point-subtracted 3.6/xm maps from Dale et al. (2009). In a few cases we 
have extended the point-subtracted mask to the outskirts of the galaxy, as the H i covers a larger 
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area than the LVL aperture in which point-subtraction was initially performed. We place the maps 
on the same astrometric grid as the H i data and convolve to our 200 pc resolution. 

To calculate the total stellar mass, -Mgtar, we sum the contribution from only the pixels 
above our Hi S/N threshold and multiply by the de- projected physical area of a single pixel. To 
calculate the global average (S*), we divide by the de-projected area covered by pixels above 
our S/N threshold. 



2.3.8. Baryonic Mass and Surface Density 

We combine our Sgas maps with the maps to find the total baryonic surface density, 
^baryon = 1.36 Shi "I" S^. We Calculate the aperture-matched baryonic mass, Mbaryon; by summing 
only the pixels above our Hi S/N threshold times the de- projected area of a single pixel. We note 



that this mass is typically smaller than the Afbaryon.tot used as a proxy for halo mass in § 2.3.2 due 
to eliminating low S/N pixels in this measurement. The global average (Sbaryon) is calculated by 
dividing Mbaryon by the de-projected area covered by pixels above our S/N threshold. 



2.3.9. Global R I Second Velocity Moment 

To facilitate comparisons with velocity dispersions in the literature, which often use the second 
moment map as a proxy, we use the intensity- weighted Hi velocity dispersion (second moment) 
maps to calculate a global second moment for the entire galaxy: 

(o-ma) = — • (6) 

id 

for every pixel above our Hi S/N threshold. We weight each pixel's second moment value, 
(Tm,2,i,i; by the Hi surface density SHi,ij. We include flux weighting in this calculation because the 
CO- added profiles discussed in § [3] are also flux- weighted. It therefore allows for a more meaningful 
comparison between the velocity dispersions derived from the superprofiles presented in this paper 
and those derived from the second moment maps as in the literature. 



2.3.10. Correlations between properties 

A number of the above physical properties are correlated with each other. As discussed further 
in § [5j we use the Spearman rank correlation coefficient, r^, to determine whether two properties 
are correlated. This statistic also yields ps, the probability of a random sample having an value 
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that is equal or more extreme than the measured value. We choose ps < 0.01 as a conservative 
threshold for correlation. 

In Figure [T| we show the correlations between many of the above physical properties. Within 
one panel, each point represents the globally- averaged properties for a single galaxy. These points 
are then colored either black if the two properties are significantly correlated or grey if they are 
uncorrelated. The correlation coefficient, r^, is shown in each panel. It is immediately clear that 
many of the mass tracers (Mhi, M^, Mbaryon, '"^20) are correlated, as expected. The total SFR is also 
strongly correlated with mass for the reason that more massive disk galaxies simply tend to have 
more material for star formation. We also see correlations among many of the more local surface 
density quantities. The correlation between (Esfr) and (Shi) is expected due to the Kennicutt- 
Schmidt relation, although this relation begins to break down in dwarf galaxies for gas masses 



measured only using Hi (e.g., Bolatto et al. 2011). The baryonic surface density (Sbaryon) also 



correlates with these quantities; given that the galaxies are gas rich with a median gas fraction of 
/gas = 0.74, (Shi) makes a large contribution to (Sbaryon)- 



3. Global Hi Superprofiles 

The global properties of the H i velocity dispersion are not necessarily well-characterized by H i 
second moment maps, as the second moment can be artificially increased by bulk motions of small 
amounts of gas at anomalous velocities. Instead, we co-add individual line-of-sight profiles after 
removal of the rotational velocity. This method produces an average, high S/N Hi line spectrum, 
which allows us to characterize the average velocity structure of the ISM. 

The basic outline of the procedure is as follows. We first measure the rotational line-of-sight 



velocity from each profile using the standard, 200 pc resolution data cube (§ 3.1). After applying 



a S/N cut (§ 3.2), we recenter each line-of-sight profile in the flux-rescaled data cube so that the 
peak is at zero and then sum all recentered line-of-sight profiles into a single, global superprofile 
(§[3l. 



3.1. Determining the Peak Velocity 

To calculate the H I superprofiles, we must first find the velocity by which to shift each line-of- 
sight spectrum in the data cube. For undisturbed, idealized Hi line-of-sight profiles, this velocity 
is simply the velocity where the spectrum reaches its maximum. However, non-circular motions 



and instrumental effects can influence the location of the peak. Initially, Braun (1997) determined 
this velocity by simply finding the velocity at which the line profile reached its maximum. In the 
comparatively higher S/N regime of Braun's data, the peak is unlikely to be strongly affected by 
noise. The peak position is also affected by the velocity resolution, such that observations with 
coarse velocity resolution cannot be used to identify the peak to better than the velocity resolution. 
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In lower S/N spectra with high velocity resolution, however, the peak of observed profiles can be 
artificially shifted either to neighboring channels or even to a completely arbitrary value by noise 
spikes. Other median line profile studies used the velocity field derived from the first moment map 
(e.g.. Dickey et al.|1990 Boulanger k, Viallefond|1992 ) or from single Gaussian fitting (Kamphuis & 
Sancisi||1993 Rownd et al.||1994 ) to determine the velocity shift. These velocity fields, however, are 
often affected by asymmetric Hi line-of-sight profiles. Because any offset in velocity can introduce 
artificial broadening into the superprofile, we must find a more robust method for determining fpcak 
by using information from the entire line-of-sight spectrum. 



Toward this goal, de Blok et al. (2008) tested a variety of methods to determine fpeak for 
rotation curve calculation: the intensity-weighted mean velocity; the velocity of the peak fiux; a 
single or a multiple Gaussian fit; and a Gauss- Hermite polynomial fit (e.g., van der Marel & Pranx 



1993). They concluded that the most robust function in the low S/N regime is a Gauss-Hermite 
polynomial that includes an /13 term. This method has already been used to generate velocity fields 



used in rotation curve analysis for a number of galaxies (e.g., Noordermeer et al. 2007 de Blok 



et al. 2008). Double Gaussian functions also tend to accurate fit line-of-sight spectra, but must 



be subjected to stringent parameter and S/N constraints to avoid fitting noise spikes (e.g., Warren 



et al. 2012). 



In individual line-of-sight spectra, we find that the first moment, single Gaussian fits, and 
the velocity of the peak can be strongly influenced by asymmetries or noise, as result similar to 



that found by de Blok et al. (2008). We show a comparison of these various peak-determining 



methods for two individual line-of-sight spectra in Figure [2} Both line-of-sight spectra show clear 
asymmetries, and the Gauss-Hermite polynomial best approximates the peak of the line-of-sight 
spectra. In comparison, the single Gaussian fit and the first moment value are both shifted due to 
the asymmetry. Because WG cirG intGrcstGci only in dctGrmining 'Upeak 

and not the detailed underlying 

structure of each individual line-of-sight spectrum, we adopt Gauss-Hermite polynomials when 
deriving Vpeak- 

We use the unmasked, standard data cube to generate fpeak maps for each galaxy. This cube 
provides the correct, uniform noise properties necessary for Gauss-Hermite fitting, even though 
the fluxes are not accurate. Nonetheless, we expect Upeak to be the same in both the standard 



cube and the flux-rescaled cube. As explained in Ott et al. (2012) and Walter et al. (2008), 



and references therein, the correction applied to the flux-rescaled cube rescales the intensity of 
the residuals of deconvolution to the same beam area as the intensity measured from the clean 
components. However, the Hi emission in channels near the peak of the profile is primarily in the 
clean components and is not rescaled. The low level of fiux in the residuals should not affect the 
location of the highest intensity emission that determines fpeak- 

To generate the Vpeak maps for each galaxy, we fit a Gauss-Hermite polynomial with an /13 
term to each pixel in the unmasked, standard data cube as given by: 



,) = Ae-y 



1 



V6 



2V2y^ - 3^/2y) 



(7) 
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where y = {v — fi)/ ctgh and the component measures an asymmetric deviation from a Gaussian 



with amplitude offset ^u, and standard deviation ctgh (van der Marel & Franx 1993). We note 
that the value of ctgh is not necessarily the same as the standard deviation of a best-fit Gaussian 
without the term. The Gauss-Hermite polynomial fits are not intended to be a representation 
of the underlying H I distribution, and we attach no physical significance to the parameters other 
than the determination of Vpeak- 

For the fitting process itself, we use a Python implementation of the Levenberg-Marquardt 
fitting algorithn|^ with uniform weight given by 1/ cJchan on each channel. We require that the peak 
velocity be within ±20 km s~^ of the first moment to ensure that the peak falls in the range of 
true H I emission. The width of the profile is forced to be greater than the velocity resolution to 
prevent the algorithm from fitting individual noise spikes. While individual line-of-sight profiles 
can show double peaks indicative of expanding structures, recent studies of spectra in these galaxies 
have shown that they comprise only a small fraction of all spectra ( Warren et al.|[2"012 ); thus, they 



do not contribute large amounts of flux to the superprofiles. We assess the accuracy of our Vpcak 



measurements in ^3.2 We further discuss the the effects of fpeak uncertainties on our results in 
Appendix |B.1[ 

As an example, the final Vpeak map for Sextans A is shown in the upper panel of Figure |3[ 
The middle and lower panels illustrate the differences between the Gauss-Hermite determination 
of 

^peak Sind what would have been determined by fitting a single Gaussian (middle panel) or using 
the first moment map (lower panel). In each map, we have overlaid a line to indicate where S/N 
> 5. The absolute differences between the Gauss-Hermite velocity field and the first moment and 
Gaussian velocity fields are small in a global sense (i.e. ~ 4 km s~^). However, the differences are 
often only a factor of two smaller than the second moment values themselves, which could lead to 
spurious broadening of the global superprofiles if we use either the first moment map or a single 
Gaussian fit to determine fpeak 



3.2. Signal-to-Noise Threshold Selection 

To ensure that we have accurately measured Upeak for each line-of-sight spectrum we include 
in the superprofile, we select only pixels for which the peak of the profile fit in Equation [7] is above 
a specific S/N threshold. For each pixel, we define S/N as the ratio between the maximum of the 
Gauss-Hermite polynomial fit and (Xchan- Figures [4] and [5] show a number of individual line profiles 
with various S/N from NGC 2366 and Sextans A, with the Gauss-Hermite fits overlaid in red. It is 
qualitatively clear that the Gauss-Hermite polynomials provide better fits as the S/N increases. For 
spectra with S/N > 5, the fitting routine does a good job at identifying the peak. For spectra with 
3 < S/N < 5, the fitting routine does a reasonable job, and spectra with S/N < 3 have decidedly 



^mpfit.py; available at 



http://code.google.eom/p/astrolibpy/ 
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questionable fits. 

To quantify tliis befiavior, we run Monte Carlo (MC) simulations for four representative galax- 
ies that span the range in velocity resolution and S/N (GR 8, Sextans A, UGC 4483, and NGC 2366). 
We create a simulated data cube from the Gauss-Hermite polynomial fits to the data, add noise at 
the appropriate level, and run our iipeak-finding algorithm. We then repeat this process 100 times. 
While the Gauss-Hermite fits are not necessarily representative of the underlying H i distribution, 
they provide a known input for our tests. The output Vpeak maps from all the MC realizations allow 
us to calculate the average uncertainty in measuring fpeak as a function of S/N. In Figure [6| we 
show the standard deviation of MC Vpcak offsets versus S/N for the four test galaxies. Each point 
represents the standard deviation around the input Vpeak of the 100 repeated Vpeak measurements 
for a single pixel. For clarity, only 5000 random points are shown for each galaxy. Regions below 
the cube's velocity resolution Av are shown in grey. 

At smaller velocity resolutions, we are better able to determine Upeak, for a fixed S/N. We 
overlay a vertical dashed line at S/N = 5, where the standard deviations of the coarsest velocity 
resolution data (2.6 km s~^) start to flatten. In all cases, we reproduce the input i^peak with an 
error of < 2 km s~^. Based on these tests, we adopt a S/N > 5 as our threshold. 



3.3. Co-addition of Line-of-Sight Spectra 

We use the flux-rescaled cubes to generate the final global superprofiles. These cubes provide 
correct flux properties and have been corrected for primary beam attenuation. The advantages of 
using the flux-rescaled cubes is that they allow us to calculate accurate estimates of H i mass and 
energy, at the cost of having more complicated noise properties than superprofiles generated from 
standard cubes. 

We apply a mask to the data cubes so that only channels with real Hi emission contribute. 



We use the 45" resolution masks described in § 2.2 for each galaxy but extend them by 15 km s ^ 



on either side in velocity. This extension includes any low-level Hi emission from gas in the 
surrounding channels that is below our masking threshold, but mostly eliminates spurious signals 
from instrumental effects such as sidelobes and clean bowls, which often occur further from the 
true Hi emission in velocity space. The final superprofiles are not strongly changed when using 
unmasked spectra in data cubes that do not show detectable instrumental artifacts such as sidelobes 
or negative bowls due to missing short spacings. 

To create the final global line profile from the masked, fiux-rescaled data cubes, we first recenter 
the selected (S/N > 5) individual line-of-sight spectra such that fpeak is at zero. Because Vpeak 
is often located in the middle of a channel, simply shifting all profiles by an integer number of 
channels and then summing can artificially broaden the final global line profile. Instead, we linearly 
interpolate across each line-of-sight profile by a factor of 10 before shifting in velocity space. Finally, 
we co-add the shifted line-of-sight profiles with equal weight to obtain the intensity-weighted Hi 



- 18 - 



superprofile. 

The final global line profiles are shown in Figure [7| ordered by decreasing Mbaryon,tot • They 



will be discussed in detail in § 4.1 



3.4. Uncertainties 

We define the noise on each point of the superprofile, asp, as: 

O'SP = 0"chan X A/^pix/A^pix/bcam X ''"'"'''''"^ . (8) 
' ^standard 

Here, o"chan is the rms noise in a single channel. A'^pix is the number of channels contributing 
to a superprofile point, and A^'pjx/bcam is the number of pixels per resolution element; this term 
represents the approximate number of independent profiles contributing to each superprofile point. 
We count each interpolated point contributing to a single superprofile point as one pixel. The 
-^rescaied /-^standard term is the flux ratio between the total measured flux in the superprofile generated 
from the rescaled cube to that from the standard cube. It is included to approximate the rescaling 
process, and maintains the same fractional noise between the standard cube and the rescaled cube. 
Typical values for FrcscaicdZ-Pstandard are 0.4 lb 0.13. We discuss the details of this estimate further 
in Appendix [A} 

The noise for each superprofile is shown in Figure [7] as the grey shaded region around each 
measured superprofile. In some cases, the S/N of the final superprofile is high enough that the 
uncertainties are smaller than the black line showing the superprofile itself. 



4. Characterizing the Superprofiles 



In this section we discuss our analysis of the superprofiles. We first give an overview of their 
general properties (§ |4.1[). We then discuss the parameterization chosen to characterize the super- 



profiles and to quantify the observed asymmetry (§ 4.2 ). We then discuss the physical interpretation 
of the parameterization (§ 4.3 ) and provide a comparison with other studies of H i kinematics (§|4.4[). 



4.1. Overview of the Final Global Line Profiles 

In Figure [7] we show the superprofile for each galaxy, overlaid with a Gaussian scaled to the 
amplitude and to the half-width half-maximum (HWHM) of the superprofiles. We also shade in 
transparent red the regions in the wings where there is more Hi than expected compared to the 
Gaussian. The superprofiles show similar structures from galaxy to galaxy — namely, a central 
narrow peak with additional wings to either side. 
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To compare the overall shape of the superprofiles from galaxy to galaxy, we plot all the super- 
profiles together in Figure |8j We normalize each superprofile so that its maximum flux is 1, and we 
scale the velocity axis by the HWHM of the superprofile. Regions where the scaled superprofiles 
overlap are darker. We also overplot a Gaussian with the same scaling, shown as the dashed line. 
Residuals from the scaled Gaussian are shown in the bottom panel. Because we have plotted each 
individual line with some transparency, regions where the superprofiles overlap are darker. 

The superprofiles in Figure [8] exhibit remarkably similar shapes, especially in the central re- 
gions. Typically, the profiles are peakier than a Gaussian in the central regions, and show wider 
wings whose residual amplitude peaks at approximately 2 x HWHM. The amplitude of the non- 
Gaussian wings varies from galaxy to galaxy, but the general shape does not change. 



4.2. Superprofile Decomposition 

We first parameterize the superprofiles with a single Gaussian. Although a Gaussian is not 
the optimal match to the detailed shape of the profile, it is a widely used parameterization and 
provides an estimate of the average profile width. Because the Gaussian shape is a poor match to 
the overall profile, especially in the wing regions, we do not perform a traditional x'^ minimization 
of the fit, but instead scale the width of the Gaussian to match the HWHM of the superprofile and 
the amplitude to match the peak of the superprofile. If we had instead fit the superprofile with 
uniform or noise-based weighting, the width of the Gaussian would increase to compensate for the 
wings. In contrast, the HWHM scaling provides a simple estimate of the average Hi kinematics 
without relying on fitting details. 

The HWHM-scaled Gaussian is shown in Figure [7] as a dashed red line overlaid on the super- 
profile for each galaxy. We measure four parameters using this HWHM-scaled fit. 

First, we measure (Tcentrai) the width of the central peak scaled to match the HWHM of the 
superprofile. This parameter characterizes the average Hi line width in each galaxy. The use of 
^central instead of the HWHM value is chosen to facilitate comparison with other studies, which 
often describe line widths in terms of a Gaussian a. We find that the median ciccntrai = 7.7 km s^^ 
with interquartile range 7.2 — 8.5 km s~^. These widths typically are > 3 times larger than the 
coarsest velocity resolution of the observations. 

Next, we measure the fraction of Hi in the wings: 

E ^S{v) - G{v)] 

f _ bl>HWHM 

/wings — ) {") 

\v\>0 

where v is the offset velocity, S{v) is the superprofile, and G{v) is the single Gaussian scaled to 
the superprofile HWHM. The /wings parameter measures the fraction of gas moving at velocities 
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faster than expected compared to the bulk of Hi. We find median values of /wings = 0.11 with an 
interquartile range 0.1 — 0.13. The measured /wings values are quite small, which indicates that the 
majority of the H i is contained in the central peak. In Figure [7| the regions of the superprofiles 
that contribute to /wings are shown as transparent red regions. 

Third, we measure the root mean squared velocity of excess flux in the wings, weighted by the 
fraction of gas moving faster expected based on the HWHM-scaled Gaussian: 



2 ^ H>HWHM 

'-'^wings • \^^) 

E ^(^) - ^(^) 

|i)|>HWHM 

The (Twings parameter measures the excess energy in the wings of the profile per unit Hi mass. It 
is also equal to the characteristic velocity of excess gas. We find median values of Swings = 21.8 
km s~^ with an interquartile range 20.2 - 25.0 km s'^. These values are typically ~ 8 times the 
coarsest velocity resolution of our data, and a factor of ~ 2 smaller than the median characteristic 
inclination-corrected rotational velocity u'2o/2. In Figure [7| Swings is shown as a solid vertical red 
line on other side of the superprofile. 

Finally, we quantify the asymmetry of the residuals in the wing regions: 



E \/iS{v) - Si-v)Y 

|i)|>HWHM 

|i)|>HWHM 



(11) 



where S{—v) is simply the mirror image of S{v) around the peak. This parameter ranges between 
in the case of complete symmetry to 1 if all the excess flux was concentrated on one side of the 
superprofile. We find a median value of a = 0.22 with an interquartile range 0.17 — 0.30. 

We have also evaluated the asymmetry of the superprofiles around the peak: 



Y^^{S{v)-S{-v)f 

Oglobal = — ■ (12) 



To first order, the superprofiles are very symmetric. We find a median global asymmetry of only 
0.05 with an interquartile range between 0.03—0.08 for the superprofiles. In all cases, approximately 
80% of the asymmetry is from the wings of the profile, while regions with velocities less than the 
HWHM are considerably more symmetric. We therefore do not include agiobai in our parameters 
as it is less sensitive to asymmetries compared to a. 
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In Table |4] we report the measured quantities from our single Gaussian fits. We list (1) the 
galaxy name; (2) cJcentrai and associated uncertainty of the HWHM Gaussian; (3) Uwings and as- 
sociated uncertainty; (4) /wings and associated uncertainty; (5) the asymmetry parameter a and 
associated uncertainty; (6) the global asymmetry Cgiobai; and (7) the number of independent res- 
olution elements comprising each superprofile, A^'bcams = -^pix, s/N>5 / -^pix/beam- We show the 
distribution of cJcentrab o-wings, /wings, and a in Figure [9j 

We discuss the determination of uncertainties on these parameters in Appendix [B) 



4.3. Physical Interpretation of the Superprofiles 



The parameterization described in § 4.2 implicitly assumes that the majority of the Hi has 
a velocity dispersion that is reasonably well-described by CTcentrai- The central peak can then be 
taken to represent the average kinematics of widespread well-thermalized H i gas across the galaxy. 
Thermal temperatures implied by the measured (Tcentrai values are ~ 4, 000 — 12, 000 K, a range 
that brackets the predicted stable ISM temperature of T ~ 7000 K ( Wolfire et al.||1995 ), but with a 
larger range. Gas in this temperature range can radiate its energy away in ~ 10^ years, a timescale 
too short to replenish the lost energy from external sources. Therefore, the value of cJcentrai niay 
be better interpreted as random turbulent velocities, which decay more slowly but still require an 
energy source to maintain over the galaxy lifetime. The deviations from a Gaussian profile can 
then potentially be explained by the fact that the central peak is likely to be the sum of Gaussians 
with a range of velocity widths due to the decline of velocity dispersion of warm H i with radius 



(Tamburro et al. 2009; Warren et al. 2012). Additionally, cold Hi with velocity dispersions < 6 



km s has been identified in some of the sample galaxies along individual lines of sight, but it 



makes up a small fraction of the Hi, with typical fractions of only < 20% (e.g.. Young et al.||2003 



Warren et al. 2012). Because cold Hi is present along some lines of sight that contribute to the 



superprofiles, its small velocity dispersion could explain why the observed superprofiles are peakier 
than a Gaussian profile. 

Gas in the wings of the superprofiles can then be interpreted as localized regions where H I is 
moving faster than expected compared to the average velocity dispersion. These anomalous motions 
likely require additional energy input to drive gas from its undisturbed state into the wings of the 
profile. The superposition of this energetic gas atop the turbulent component of the central peak 
produces superprofile wings with amplitudes higher than expected from the turbulent component 
alone. 

An alternative physically-motivated decomposition of the superprofiles is to consider them to 
be the sum of the CNM and WNM, characterized by narrow and broad velocity components. This 



approach has been recently pursued in a similar study of THINGS galaxies (12012), motivated by 



the fact that H i can exist at two stable temperatures and that some individual line-of-sight spectra 
show evidence of these two phases. In this scenario, the ratio of H i flux in the narrow component to 
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that in the broad component provides a measurement of the relative amounts of cold and warm H i. 
However, this interpretation also implies that the two Gaussian components are well-thermalized 
but independent, and that only two distinct warm and cold gas populations exist in a single galaxy. 
It also presumes that the line widths are directly connected to the thermal temperatures of the gas, 
which is difficult to reconcile with both the short thermal timescales and the mismatch between 
the inferred kinetic temperatures and the predicted thermal temperatures of the CNM and WNM. 
We have found that double Gaussian fits to the superprofiles are indeed a good representation of 
the overall shape, but we believe that interpreting the Gaussian components as representative of 
the CNM and WNM is not necessarily convincing for the global profiles, although it can be valid 
along individual lines of sight. We discuss the reasons we have chosen not to use this method in 
more detail in Appendix [Cl 



To facilitate comparison with 12012, however, we have also fit the superprofiles with a double 
Gaussian using the 12012 methodology. As in 12012 the double Gaussian fits are weighted by the 
inverse of the approximate uncertainty due to noise on each point. The results of these fits are 
given m Table [5} We list (1) the galaxy name; (2) the width of the narrow component, (t„; (3) the 
width of the broad component, at', (4) the area of the narrow component relative to that of the 
broad component, An/Ab; and (5) the width of the narrow component relative to that of the broad 
component, (Tn/cTb- 



For the reasons given in § 4.2, 4.3, and Appendix O we have chosen to parameterize the 



superprofiles based on the simple HWHM scaling discussed above. We interpret the superprofiles 
physically as exhibiting a central turbulent peak with more energetic gas in the high- velocity wings 
to either side. 



4.4. Comparison with Other Studies 



The H I velocity dispersion and associated energy have traditionally been estimated either by 
fitting single Gaussian profiles to line-of-sight spectra (e.g., Petric & Rupen 2007) or by using 



the second moment map (e.g., Tamburro et al. 2009). Since our method is an uncommonly-used 



estimate of H I turbulent velocity, we provide comparisons between our work and previous studies 



from Petric Rupen (2007) and Tamburro et al. (2009) 



Petric & Rupen ( 2007 ) determined cjhi by fitting single Gaussians to each line-of-sight spectrum 
in the face-on spiral NGC 1058 at a spatial resolution of ~ 1.3 kpc, much larger than the 200 pc 
scales studied in this paper. The typical range of velocity dispersion in the disk was 4-14 km s~^, 
with a majority at ~ 7 km s^^. They also found that the median profile shape, after normalization 
for line-of-sight profile width, was similar to the shape observed in our superprofiles, characterized 
by a central narrow peak with wider wings. Many other studies have found similar shapes for the 
average Hi line profile (e.g.. Dickey et al.||1990 Boulanger Sz Viallefond|[l992 Kamphuis & Sancisi 



1993), but quantitative comparisons with these studies are hampered by large differences in velocity 
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resolution or physical resolution. 



Tamburro et al. (2009) used the second moment as a a proxy for Hi velocity dispersion. They 
found that galaxies have turbulent components with amplitudes of ~ 10 km s^^. When we compare 
the four galaxies that overlap both samples (Ho H, NGC 4214, IC 2574, NGC 7793), we measure 
consistently lower ^central values than the global second moment by ~ 2 km s~^. We do not believe 



that the difference results from the fact that the Tamburro et al. (2009) sample focused mainly 
on large spirals compared to the low- mass dwarfs that dominate our sample, as all galaxies in our 
sample have cJcentrai values that are smaller than their global second moments. The difference is 
more likely due to the fact that we have isolated the central peak in our sample, while the second 
moment values are affected by the presence of high velocity wings and asymmetries in the line 
profiles. 

In Figure 10 we show a comparison between our measured cicentrai value and ((Tm2) (Equation[6|, 
which provides an estimate of the turbulent width that would have been derived using the |Tamburro| 



et al. (2009) methodology. Compared to cTcentrah iii all cases the global second moment is larger 
by ~ 10 — 50%, with a median of 20%. When interpreting the line widths physically, this offset 
suggests that the second moment leads to higher estimates of the energy necessary to drive Hi 
turbulence. 



We can also compare our results to those of 12012, who have performed the most directly 
analogous analysis to date. In this study, the authors generated superprofiles for a number of 
THINGS galaxies, eight of which overlap with our sample. However, their approach was somewhat 
different. The authors used the naturally-weighted standard maps and worked at the instrumental 
resolution. They also chose to model each superprofile as the sum of a narrow and a wide Gaussian 
profile representing the CNM and the WNM. 

The first difference to note is that our sample is smoothed to a fixed physical resolution, as 
opposed to working at the instrumental resolution of ~ 10" for each galaxy. While smoothing to a 
coarser resolution can broaden the intrinsic superprofile because the measured velocity dispersion 
increases at larger physical scales, matched-resolution cubes allow a more robust comparison from 
galaxy to galaxy by sampling the Hi kinematics on the same physical scale. Since the more massive 
THINGS galaxies tend to be at larger distances, the varying spatial resolution could lead to spurious 
trends with any quantity that correlates with galaxy mass. 



The superprofiles that we derive are systematically different from 12012 even if we fit them 
with the double Gaussian function. For the nine THINGS galaxies in common (DDO 53, DDO 
154, Ho I, Ho H, IC 2574, M81 DwB, NGC 2366, NGC 4214, NGC 7793), we typically measure 
similar (T„ values, but find systematically smaller values of at, by 0.5 — 3 km s~^. This offset is also 
apparent in the (Tn/cTb ratio; we find an average OnjcTh = 0.46, nearly 25% higher than measured 
bypOT2 



We believe that this difference could result from generating our superprofiles from the flux- 
rescaled cubes instead of the standard cubes; from using robust-weighted cubes instead of natural- 
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weighted cubes; or from choosing different S/N thresholds. 

First, the use of flux-rescaled cubes hkely lowers the amount of flux in the wings of the profile. 
While the central peak of the superprofile is due to the regions of the line-of-sight spectra that 
are the brightest, the wings arise from low-level flux to either side of the peak. When we consider 
an individual line-of-sight Hi spectrum, the peak of the spectrum has a higher fraction of flux in 
the clean components than in the residuals when compared to the lower-level emission to either 
side of the peak. Because the flux-rescaling correction effectively lowers the amount of flux in the 
residuals, the low-level emission on either side of the peak in each line-of-sight spectrum is rescaled 
to a relatively smaller value compared to the peak regions. The flux-rescaling correction therefore 
has the effect of narrowing line-of-sight profiles, which then produces a narrower central peak in 
the superprofiles and less fiux in the wings compared to superprofiles generated from the standard 
(non-rescaled) data cube. With smaller wing amplitudes, the broad Gaussian component a^, is 
smaller, which could then make the Onloh ratio larger. 



Second, the synthesized beam of natural-weighted data cubes, which were used by 12012 
exhibits a positive pedestal that extends to large radii which is not present in the robust-weighted 
data cubes that we have used to generate the superprofiles. Because of the broad pedestals in 
the synthesized, naturally-weighted beams, each apparently independent beam, as judged by its 
FWHM, actually includes flux from a much greater area. These pedestals can therefore lead to 
additional velocity smearing, beyond what one would expect based on the velocity fleld of the 
galaxy and the angular size of the beam. The neighboring pixels that are included at low-level, 
however, have offset velocities due to the overall rotation of the galaxy. Therefore, each line-of-sight 
spectrum also includes flux from neighboring pixels at offset velocities. This additional fiux at offset 
velocities may translate to more fiux in the wings of the superprofile, and therefore larger measured 
values of o\, and smaller (Jnlcyh ratios when compared to the robust- weighted cubes that lack these 
positive pedestals. 

Finally, the signal-to- noise threshold was different between the two studies. At lower S/N 
values, determination of the peak is more difficult (see Figure [g]). Because the determination of 
I'peak is worse for low S/N spectra, additional fiux could be incorrectly added into the wings of the 
superprofile due to these offset Vpeak values. This addition of fiux into the wings would likely widen 
the cJb measurement and lower the On / value. 



5. Comparison with Physical Properties 

In this section we examine how global physical properties of the galaxies correlate with the 



properties of the global Hi superprofiles, as characterized by (Jcentrah cr„ings, /wings, and a (§ 4.2). 



For all tests, we characterize the strength of the correlation compared to a random sample 
using the Spearman rank correlation coefficient. This method makes no assumptions about the 
functional relationship between the two input data sets, and instead tests only for a monotonic 
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relationship between the two variables. The Spearman correlation coefficient varies between 
-1 (monotonically decreasing) and 1 (monotonically increasing) with implying no relationship. 
The significance of the observed value of is given by where ps < 0.01 means that random 
uncorrelated data produce a value at least as extreme as measured < 1% of the time. We adopt 
this threshold to indicate a statistically significant relationship between two quantities. 

For all following comparisons, we list and ps in Table [6j The coefficients for significant 
correlations [ps < 0.01) are shown in bold. 

Throughout this section we also present a number of figures showing the behavior of the 



superprofiles with a different physical property, starting with Figure 11 The four upper panels 
show cTcentrah (Swings, /wings, and o. Colors are determined by the physical property itself, with 
red indicating a low value and blue indicating a high value. In the lower panel, we show the 
scaled superprofile residuals from the HWHM Gaussian fit, as previously seen in the lower panel of 
Figure [8| In this case, though, we have added color to each line to highlight how the superprofile 
shape changes with that physical property. The colors of each line are the same as the corresponding 
point in the upper panels, though the lines have a transparency value so that overlapping regions 
are clearer. To better highlight any asymmetry of the superprofiles, we have reversed the velocity 
axis if necessary such that the wing peak with the higher area is on the left. We now discuss these 



figures for correlations with the physical properties calculated in § 2.3 We also remind the reader 
that many of these properties are physically connected with each other, as shown in Figure [T| so a 
correlation with one property could be causally due to another. 



5.1. Correlations vi^ith GalELxy Mass and Related Quantities 

We start by examining trends in the superprofile parameters versus quantities that correlate 
with galaxy mass. While we would not expect the local Hi conditions to know much about the 
overall galaxy potential, galaxy mass is correlated with a host of other properties that are more 
likely to have direct effects on the Hi, such as SFR or Mhi. 



5.1.1. Halo Mass 



As a proxy for galaxy halo mass, we use M\y, 



iaryon,tot 



(§ 2.3.2). In Figure 11 we plot the 



correlations between galaxy mass cind. the parameters (Tcentrab 'Swings ^ /wings, and a derived from 
the superprofiles. Neither ^central nor o"wings shows any significant trend with galaxy mass. This 
lack of correlation confirms our expectation that the driver of H i kinematics is not strongly mass- 
dependent; that is, it is more likely to be dependent on specific phenomena in the ISM. 



However, Figure 11 does show a strong correlation between galaxy mass and the asymmetry a, 
with more massive galaxies exhibiting less asymmetry. The asymmetry can be seen in the residuals; 
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the red lines of low-mass galaxies show a large variation between one side of the residuals and the 
other, while the blue lines of higher-mass galaxies have more symmetric structure in the residuals. 
However, it is likely that this trend reflects of galaxy properties such as SFR that correlate with 
mass, instead of mass itself. Asymmetries in the line-of-sight spectra induced by star formation 



may propagate to asymmetries in the superprofiles, as discussed further in § 6.2.3 but are averaged 
out when more star forming regions are present. 

We also find a trend between galaxy mass and /wings- The ps value is 0.017, only slightly above 
our cutoff for a significant correlation. There is a tentative indication that lower mass galaxies with 
-^baryon.tot ^ 5 X 10® Mq tend to have a larger scatter in /wings- Galaxies with larger Mbaryon,tot 
exhibit higher /wings values with less scatter, showing half the standard deviation relative to the 
median compared with their lower mass counterparts. If this trend holds true with a larger sample, 
it would be consistent with higher mass galaxies being able to more consistently drive H i into the 
wings of the superprofile. The presence of high velocity gas in lower mass galaxies could be due 
to more stochastic processes, with not every galaxy being able to launch high velocity gas at all 
times, thus leading to more scatter in /wings - 

We have also looked for correlations with other mass tracers, such as Mhi, M^, and W2q- While 
the exact values of and ps change, we find similar trends between superprofile properties and 
other mass tracers. Such behavior is expected expected based on the strong correlations between 
Mbaryon.tot ^-hcI other mass tracers shown in Figure [T] 



5.1.2. Star Formation Rate 

We plot the behavior of the superprofile parameters with SFR in Figure [12] As shown in 
Figure [T| SFR is strongly correlated with galaxy mass. We are therefore not surprised to see 



nearly-identical correlations as those in Figure 11 We do find statistically significant correlations 
between SFR and both /wings and a. The trends can be seen in the upper panels of Figure 12 where 
galaxies with higher SFRs (blue) tend to have higher /wings values and lower o values compared to 
the galaxies with lower SFRs (red). This behavior can also be seen in the superprofile residuals 
plotted in the lower panel of Figure [T2| the blue residuals, with the highest SFR, are also among 
the largest and most symmetric. The red lines are typically lower and exhibit varying levels of 
asymmetry. 

Given that the correlations between /wings and o are the same as the correlations with mass, 
it is possible that the /wings and a correlations are causally connected to any of the other galaxy 
properties that scale with mass. However, of all these properties, SFR is the only one that provides 
a physical mechanism for driving Hi gas in the wings. It may therefore be the actual driver of the 



correlations with mass. We explore this connection further in § 5.2 below. 



-27- 



5.2. Star Formation 



Star formation is typically proposed as the primary driver of H i turbulence in spiral galaxies 
within ~ (e.g., Tamburro et al. 2009). If so, then we may expect to find correlations between 



measures of star formation and the superprofile parameters. While the total SFR provides a 
measure of the overall energy input from star formation, it is strongly dependent on galaxy mass. 
In addition, galaxies with higher SFRs, and thus higher masses, also have larger H I masses to affect 
with the energy provided by star formation. 

We consider two possible measures of star formation other than SFR, which was considered 
previously in § 5.1.2, These include the SFR intensity (i.e., (Ssfr); § 2.3.5), which measures the 



SFR concentration and thus may correlate with the efficiency of locally accelerating H i; and SFR 
/ Mhi (§ 2.3.6), which measures the ratio between the available energy of star formation and the 



mass of gas that the energy must couple to. 



5.2.1. SsFR 

In Figure 13 we show the relationship between (Ssfr) and the superprofile parameters. We 



find no significant correlations with o"ccntrah /wings, or a. However, fiwings shows a trend with (Ssfr), 
such that galaxies with higher (Ssfr) values are able to drive Hi into the wings with faster average 
velocities. 

The existence of the correlation between Uwings and (Ssfr) is not surprising. Since energy input 
from star formation is a local process, more concentrated star formation should be more effective at 
inducing anomalous motions in the surrounding Hi. At the same efficiency of converting the star 
formation energy into kinetic H I energy, a higher concentration of star formation energy imparts a 
given amount of energy into a smaller mass of H I, thus driving H I in the wings to higher velocities. 
We further explore the connection between star formation energy and kinetic energy in the wings 



of the superprofile in § 6.3.4 



5.2.2. SFR / Mhi 



We might also expect SFR per unit Mhi to affect Hi kinematics. As discussed in § |2.3.6[ this 
quantity is best interpreted as the ratio between the rate of energy input from star formation and 
the mass of gas that can be accelerated by that energy instead of the average SFE. 



We show the superprofile parameters as a function of SFR / Mhi in Figure 14 We find a 
correlation between /wings and SFR / Mhi. As a potential explanation for this trend, we invoke 
the argument that galaxies with higher SFR and lower Mhi are better able to accelerate the 
surrounding H I into the wings of the superprofile compared with their counterparts. However, for 
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the most common values of SFR / Mhi ~ 10 yr ^, /wings values vary by ~ 40%, suggesting a 
large degree of stochasticity in this correlation. It is also unclear why SFR / Mhi affects /wings but 

not Cwings- 



5.3. Surface Mass Density 



Disks with higher surface mass density are able to permit more turbulent motions in gas that is 
still bound to the disk (e.g., van der Kruit|198T ). In addition, a number of gravitational instabilities 
have been proposed to drive turbulence in galaxy disks (e.g., Huber &: Pfenniger||2001 Wada et al. 



2002 ) , so the amplitude of turbulent motions might be expected to scale with surface mass density. 



We find correlations between (Shi) and both (Jcentrai and CTwings, as shown in Figure 15 This 



correlation may be indicative of turbulence as driven by gravitational instabilities, as explored 
further in § 6.3.1, or of a coupling with (Ssfr); which tends to scale with Shi. We note, however. 



that no correlation was observed between cJccntrai and Ssfr- A similar correlation with (Sbaryo 

exists with CXwings, but not with Ucentral- 



5.4. Inclination 

We end by confirming that the superprofile parameters are not strongly influenced by galaxy 



inclination, as can be seen in Figure 16 There are no significant trends in any of the parameters. 
However, we note that uncertainties in the inclination of the dwarfs could potentially mask any 
underlying systematic effects. 



5.5. Extending the Correlations to Higher Mass Galaxies 



A number of massive spirals were eliminated from our sample based on the selection criteria 



to ensure high-quality data (§ 2.1.3). In this section we now include eight of these galaxies to see if 
the correlations we identify among the dwarfs could hold when extended to higher mass galaxies. 
Since this is a simple check, we do not perform the same rigorous tests as we have for the primary 
dwarf sample. Instead, we merely present our results and assess how the correlation coefficients 
change with the inclusion of more massive galaxies. 



For the higher mass galaxies, we relax selection criteria (1) and (2) as given in § 2.1.3 by working 
at a spatial resolution of 400 pc and by including observations with At; = 5.2 km s~^, with the 
caveat that we have not characterized the effect of observational properties on these superprofiles. 
We also include galaxies with de Vaucouleurs T-type > 2. Finally, we only use galaxies that are in 



the clean sample of 12012 to avoid contamination from bulk inflows, outflows, or Hi in the Milky 
Way. These relaxed criteria give us eight additional, higher-mass galaxies (NGC 628, NGC 2403, 
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NGC 2903, NGC 2976, NGC 3351, NGC 4736, NGC 5055, and NGC 5236). 

For each galaxy, we produce superprofiles in the same manner as we have for the dwarfs in 



our sample. We first smooth the data to a common physical resolution of 400 pc (following § 2.2). 
We then generate a superprofile for that galaxy after applying the same S/N > 5 threshold as 
for the primary sample (§ [s]). Finally, we parameterize each superprofile using the HWHM-scaled 
Gaussians (§ 4.2). We list the derived superprofile parameters for these galaxies in Table [Tj 



We assess the strength of the correlations with physical properties after including the higher- 



mass galaxies. In Figure 17 we show relevant correlations between the superprofile parameters and 
the physical properties with which they are correlated after the inclusion of high-mass galaxies. The 
displayed panels are: (1) iTcentrai versus (Shi); (2) cJcentrai versus (Ssfr); (3) cr„ings versus (Shi); 
(4) cTwings versus (Ssfr); (5) SFR / Mhi versus /wings; and (6) SFR versus a. The dwarf galaxies 
from our sample are shown with filled circles, while the higher mass galaxies are shown with open 
circles. 

We find that o"centrai is no longer correlated with (Shi), but a trend with (Ssfr) is present 
with rg = 0.404 and ps = 0.024. Although this ps value is marginally higher than our cutoff, the 
trend lends credence to the idea that the observed correlation between (Tcentrai and (Shi) for the 
dwarf sample may be tracing a correlation with Ssfr. 

Similarly, Uwings no longer shows correlations with (Shi) but does with (Ssfr), implying that 
Ssfr does indeed affect the gas in the wings. The idea that H l gas in the wings of the superprofile 
is influenced by star formation is also supported because correlations between /wings and SFR / 
Mhi and between a and SFR both remain with the inclusion of higher mass spirals. 



6. Discussion 



In this section we discuss the implications for the width of the central peak of the superprofile 
(§ 6.1) as well as the correlations we see with the wings of the superprofile (§ 6.2). We then 



compare the energy available from star formation to the kinetic energy in the Hi gas (§ 
Finally, we approximate implied Hi scale heights for the sample 



6.3) 



6.4) and discuss the similarity 



of the superprofiles' shapes (§ 6.5) 



6.1. Hi in the Central Peak 



The width of the central peak, CTcentrah has a very small range in our sample, with a median of 
7.7 km s^^ and interquartile range of only ±1 km s~^. The small range of observed Ccentrai values 
may suggest that average turbulence in the ISM is regulated in some way, such as by energy input 
from physical processes or external heating from the UV background (e.g., Schaye|[2004 Tamburro 
et al.||2009D. 
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The only correlation we have found between cxcentrai and globally-averaged physical properties 
is with (Shi), but this correlation disappears with the inclusion of higher mass spirals. Interestingly, 
we find no trend between cTcentral and our measurements of star formation. In standard lore, star 
formation is the primary driver of turbulence in the warm ISM. If star formation were the sole 
driver of Hi line widths in galaxies, we would initially have expected to have seen a connection 
between some measure of SFR, (Ssfr), or SFR / Mhi. We now give some potential explanations 
for this mismatch. 

Because regions with higher (Sm) also tend to have higher (Ssfr), the correlation between 
^central and (Sri) could in actuality a correlation between (Tcentrai and (Ssfr)- However, this 
interpretation is called into question by the lack of correlation between cXcentrai and (Ssfr) as traced 
by FUV+24/im emission. The FUV+24/im tracer probes SFR averaged over the past 10 — 100 Myr, 
with the implicit assumption that the SFR has been constant over that timescale. It is unlikely 
that this is the case in our sample galaxies. Because the timescales over which turbulent gas in the 
central peak can dissipate energy are ~ 10 Myr, the FUV+24/xm timescale may be a poor match 
to the timescales relevant for the Hi component. If we interpret the central peak as turbulent, the 
Hi in the sample galaxies is able to dissipate its energy in ~ lO"^ yr. Therefore, the star formation 
that has occurred in the last ~ 10^ yr is the primary influence on the gas. Because our SFR 
measurement has been averaged over a longer time, galaxies with similar average SFRs on ~ 10^ 
yr timescales may in fact may have different SFRs in the past ~ 10'' yr. 

We caution that the superprofiles for the galaxies are a global average of individual Hi line- 
of-sight spectra, so they often include regions in a single galaxy with very different star formation 
properties. The loss of spatial information in the superprofiles could account for the lack of correla- 
tion between star formation and cJcentrai, especially when regions with prominent star formation are 
mixed with those that lack strong star formation. We explore the spatial dependence of velocity 
dispersion in a subsequent paper. 

It is certainly possible that the correlation with (Sm), and not with (Ssfr), indicates that 
there is indeed no physical connection between the central peak and star formation. There are hints 



of this possibility in the Tamburro et al. (2009) results. In the inner regions of their more massive 
galaxies, they found a correlation between a"ni2 and Ssfr- However, this correlation broke down at 
large radii, where (Tni2 approaches a nearly constant value but star formation falls off dramatically. 
Given the similarity between dwarf galaxies and the outer disks of spirals (in terms of Ssfr, Shi, 
S^, etc.), it is possible that our sample lies primarily in this regime. 

If star formation is not setting the velocity dispersion, what else is? One commonly-adopted 
mechanism in spirals is the MRI, which works well in the outskirts of massive spiral galaxies 
where angular velocity declines with radius. In dwarfs, however, the rotation curves across much 
of the observable disk are closer to solid body rotation, and therefore lack the strong differential 
rotation necessary for the MRI. Most dwarfs in the Local Group also show magnetic field strengths 



approximately three times smaller than that of spirals (Chyzy et al. 2011). The combination of 
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these two factors means that the MRI should be less effective in the dwarf galaxies that comprise 
our sample. Since more massive galaxies host conditions that are more conducive to MRI-driven 
turbulence, we might expect to see a correlation between galaxy mass and cicentrai- However, no 
such trend is present. While it is possible that the range in cicentrai values is too small to measure 
such a trend, it would also be surprising if the MRI conspired to produce such similar (Tcentrai values 
across the sample without any external regulation. 

Other energy sources for turbulence in galaxies include gravitational instabilities. Many of 
these, however, require shear in the rotating gas to function and thus fall prey to the same problems 



as the MRI. The most promising of these instabilities is that presented in Wada et al. (2002), which 



does not explicitly require shear to drive turbulent velocity dispersions and hearkens back to the 
observed correlation between cXcentrai and (Shi). While the authors have shown that this method 
can drive turbulence at levels observed in NGC 2915, others have noted that this instability provides 
energy at levels that are two orders of magnitude smaller than that required to drive the observed 



turbulence in the ISM (MacLow 2004). We assess the ability of this gravitational instability to 
provide enough energy to drive turbulence in our sample galaxies in § |6.3.2 



Another proposed method is UV heating, which can drive thermally-broadened line widths to ~ 
6 km ( Tamburro et al.|2009 ). In this case, the widths are due to thermal effects, not turbulence. 
However, the measured fJcentrai values of the superprofiles show a much wider temperature range 
than can be explained by UV heating. Perhaps, however, UV heating can sets the base velocity 
dispersion of Hi profiles, and any additional dispersion is driven by other physical processes, such 
as star formation or instabilities. 

The mechanism that drives Hi velocity dispersions remains an open question and will likely 
become clearer in future spatially-resolved studies. 



6.2. The Energetic Hi in the Wings 

Compared with cJcentrai, the properties of the superprofile wings are more correlated with galaxy 
physical properties. We found that: (1) the characteristic velocity of the wings, cJwings, increases 
with both SsFR and Shi; (2) the fraction of the gas in the wings, /wings, increases most strongly 
with SFR / A/hi; and (3) the asymmetry, a, is primarily in the wings and decreases with increasing 
SFR. In this section we discuss potential physical explanations for these trends. 



6.2.1. SsFR, Shi, and a^-mgs 

We find correlations between cxwings and both (Ssfr) and (Shi). The fact that both correlations 
exist may be due to the (Ssfr) - (Shi) connection. The correlation between (Things and Shi may 
also partially be due to the correlation between cxcentrai and Shi; galaxies with wider central peaks. 
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and therefore higher values of Shi, wiU by definition also have higher a^ngs values, as we consider 
only gas moving faster than expected compared to the central Gaussian when calculating cTwings- 

We 

therefore discuss o"wings properties in relation to star formation, a potential driver for high-velocity 
gas in galaxies. 

The behavior of Swings and Ssfr can be explained if the energy from star formation pushes 
Hi to higher velocities than expected from the ^central Gaussian. Expanding Hi structures have 



been observed in numerous studies down to the instrumental resolution (e.g., Brinks & Bajaja 
1986t [Bagetakos et al.||201lD , including many smaller structures in the Milky Way (e.g., |Ehlerova fc 
PalouspOOS ) . Presumably smaller expanding H i structures exist below the current limits of spatial 



resolution. Both Type 2 and Type 3 holes, as defined by Brinks Sz Bajaja (1986), show Hi offsets 



in velocity space along a single line-of-sight spectrum, which would contribute Hi emission to the 
wings of a global superprofile. 



Hi holes are often thought to be due to star formation (e.g., McCray Sz Kafatos|1987 ), although 



direct spatial correlation with young, massive stars is not always seen (e.g., Rhode et al. 1999). 
Recent studies have found that multiple star formation events over the age of the hole do provide 
enough energy to drive H I hole formation, though other regions show similar star formation histories 



without the presence of Hi holes (e.g., Weisz et al. 2009 Warren et al. 2011 Cannon et al. 2011). 
H I gas with anomalous velocities has also been linked to the presence of H i holes and star forming 



regions in the spiral NGC 6946 (Boomsma et al. 2008). However, projection effects make these 



measurements difficult, so such studies are only appropriate in a small number of face-on spirals. 
If similar expanding Hi structures exist at smaller spatial scales, we can explore their connec- 



tion with superprofile parameters in more detail using the canonical Chevalier equation ( Chevalier 



1974), which relates star formation energy and Hi hole properties: 



E (erg) = 5.3 x 10^^ 



no 



cm 



^'holc 



3.12 



'.'cxp 



1 km s 



1.4 



(13) 



Here, E is the enclosed, single-burst energy required to drive the expansion, uq is the density of 
the surrounding ambient medium, rhoie is the radius of the hole, and Vexp is the expansion velocity. 
This equation provides a description of the energy of an expanding shell in an idealized ISM, with 
the assumption of uniform density which does not hold true on global scales. This equation is 
also complicated by the fact that multiple bursts can provide energy to drive Hi hole expansion. 
However, it still allows us to explore the relation between the physical quantities involved in the 
expansion. The energy density of the hole, ~ -Ehoie/?^^, scales most strongly with Vexp, with weaker 
dependencies on no and rhoie- If we assume that the energy source for the hole is due to star 
formation, a higher concentration of star formation energy within the hole should lead to faster 
expansion velocities. Because there is no reason not to expect smaller expanding structures below 
the 200 pc resolution of our data, these structures may manifest as the high- velocity wings observed 
in the superprofiles. Thus, for a given Hi mass, higher concentrations of energy due to larger SgpR 
should lead to faster expansion velocities and therefore to higher measured ciwings values. It is 
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unclear, however, why Ssfr would influence cr^mgs without also driving more gas into the wings of 
the profile, as measured by /wings- 



6.2.2. SFR / Mm and Uings 

The fraction of gas in the wings is correlated not with Ssfr but with SFR / Mhi for dwarf 
galaxies. Such behavior can again be ascribed to expanding Hi structures. At a fixed SFR, and 
therefore a fixed energy input into the ISM, galaxies with smaller Hi masses should show more 
pronounced effects on the Hi kinematics. Therefore, galaxies with relatively high SFRs compared 
to their H I masses should be able to perturb the H i content more easily, as seen in the correlation 
between SFR / Mhi and /wings- More Hi is pushed into the wings of the superprofile if the galaxy 
has a high SFR relative to its Hi content, or if it has a smaller amount of Hi to move around with 
the energy available from star formation. 

The correlations between /wings and dwings indicate that star formation does indeed play a role 
in driving Hi to anomalous velocities seen in the wings of the superprofile. It is unclear why /wings 
and not Swings would scale with SFR / Mhi. However, the correlation coefficient between SFR / 
Mhi and Swings implies ps = 0.066, which indicates only a 7% probability of finding a correlation 
this extreme from a random sample. 



6.2.3. Star Formation as a Driver of Asymmetry 

Since the properties of the superprofile wings appear to be connected to star formation, it is 
not surprising that asymmetries in the wings can also be attributed to star formation. Due to 
the inhomogeneity of the ISM, individual star forming regions can affect the local Hi gas asym- 
metrically. This effect can then result in asymmetric line profiles near star forming regions, which 
contributes to asymmetry in the wings of the global superprofile. The average of a large number of 
asymmetric H I regions should average out to produce a symmetric superprofile, while the average 
of only a few asymmetric regions is more likely to retain net asymmetry in the superprofile. If 
galaxies with larger SFR have more individual star forming events compared to their counterparts 
with smaller SFR, we would expect an anti-correlation between SFR and a, as is observed. We 
note that it may also be easier for more massive galaxies to remove the signatures of asymmetric 
Hi motions, regardless of their origin, due to their deeper gravitational potential wells. 



Asymmetric Hi motions have already been observed near star forming regions. Young et al 



(2003) found that Hi line-of-sight spectra exhibit asymmetry near regions of star formation on 



200 pc scales, though this behavior was not seen in the sample observed by Begum et al. (2008), 
with somewhat more coarse resolutions of 300 - 700 pc. We also find that some galaxies show 
asymmetric profiles near regions of star formation. To measure this, we use the difference between 
Wpeak from Gauss-Hermite fits and the intensity-weighted mean velocity from the first moment map 
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(^^iwm)- As previously seen in Figure [2| viwM and fpeak from Gauss- Hermite fits are offset for 



asymmetric profiles. In Figure 18, we plot Ssfr compared to the difference between I'peak and 
viwm/ Red indicates regions where Wpeak < "^iwivi, while blue shows regions where Vpeak > "i^iWM 
km s~^. Transparency has been added to show the underlying Ssfr in grey, and regions where 
^'peak = '^IWM are more transparent than regions where they are different. In these galaxies, some 
star formation regions also show asymmetric profiles. The spatial overlap between star formation 
and asymmetric profiles is not proof that star formation is the cause of asymmetry, as there are 
regions with no apparent star formation that also show asymmetric profiles, but it indicates that 
star formation may be one cause of asymmetry in H i line-of-sight spectra. If this is the case, the 
correlation between o and SFR is not surprising. 

We next estimate the number of star forming events, A'^sf, in a galaxy based on its SFR plus 
a fiducial star formation timescale and mass. We first assume that star formation is linearly pro- 
portional to II2 mass, with a timescale rjep ~ Shj/Ssfr. This timescale estimates approximately 
how long it will take a galaxy to use up its entire II2 reservoir, and has been found to remarkably 
independent to environment ( Bigiel et al.||20lT Bolatto et al. 2011). We can then estimate Nsf- 

SFR X Tdep 



A^SF 



MsF 



(14) 



where Msf is the typical mass of a star forming region. The numerator is an estimate of H2 mass 
in our galaxies, while the denominator is the average mass of a star formation clump. 

It is becoming clear that Tdep may be universal as it does not appear to vary much from galaxy 
to galaxy (e.g., Bigiel et al. 2011). Bolatto et al. (2011) measure Tdep ~ 1.6 Gyr in the Small 
Magellanic Cloud for the same spatial scales of 200 pc as our data, a value very similar to the 



Tdep ~ 2.35 Gyr measured by Bigiel et al. (2011). The SMC value is likely the best comparison. 



since it has an Hi mass and SFR similar to our sample (Stanimirovic et al. 111999 [Harris Sz Zaritsky 



2004). If we assume that star formation arises from Giant Molecular Clouds (GMCs) with average 
masses of Mmc > whose sizes are well-matched to our 200 pc scale and therefore provide an estimate 
of the number of 200 pc resolution elements with star formation, we can calculate the number of 
spatially-resolved elements with star formation: 

Ns. = .M ^f.*! J (r^^)(.-^y- (15) 



1 X 10-3 M, 



yr" 



1.6 X 109 yr7 VI x 10^ 



This equation provides an approximation of the number of resolution elements in our sample that 
have star formation. 

At the low end of our sample, the observed SFR ~ 1 x 10"^ Mq yr"^ kpc^^ yields only a few 
regions of active star formation. Indeed, many of the low-mass galaxies in our sample show only a 
few clumps of star formation as traced by FUV + 24/im emission, while the higher mass galaxies 
have a more widespread, smooth star formation distribution across their disks. With our above 
assumption that each SF event has a chance to drive asymmetric H i outflows, we would therefore 
expect that these potential asymmetries do not always average out at the low star formation 
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rates characteristic of our sample. At the higher end of our sample, where SFR ~ 1 x 10~^, we 
would expect a few hundred regions of active star formation. While each individual region may 
produce asymmetric H i motions, the average over the large number of SF events in the entire galaxy 
produces an overall symmetric distribution. However, very large Hi holes or extreme star formation 
events could still produce observable asymmetric Hi outflows even in galaxies with relatively high 
star formation rates. This behavior has been seen in NGC 2366 (van Eymeren et al. 2009), which 
also has a high a = 0.22 value compared to its relatively high SFR, though it is unclear if this 
asymmetry is due to star formation or the high degree of non-circular motions in the northwestern 
region (lOh et al.||2011|). 



6.3. Comparison of Energy in the Hi gas to Energy Sources 



In this section we estimate the kinetic energy of the H i gas using the superprofile and compare 
it to the energy available from physical processes. In § 6.3.1 we assess the ability of both the 



Wada et al. (2002) gravitational instability (§|6.3.2) and star formation (§ 6.3.3) to provide enough 



energy to drive turbulence at levels indicated by the central peak. We then examine the efhciencies 



required for star formation to move Hi into the wings of the superprofiles (§ 6.3.4). Finally. 



we discuss whether star formation can drive the full Hi kinematics measured by the superprofile 
(§[6X5|). 



6.3.1. Energy in the central H i peak 

We first estimate the kinetic energy contained in the central peak of the superprofile. We 
assume that the majority of gas in the central peak can be reasonably approximated by a Gaussian 
profile with a width (Jcentrai and that the velocity dispersion is isotropic in three dimensions. The 
energy in the central peak is therefore: 



3 

-E'ni.central ~ ~ Swings) (1 - /cold) A^HI 0"ccntral (16) 



where Mhi is the measured Hi mass of the superprofile. The (1 — /wings) factor accounts for the 
fraction of Hi in the wings, so (1 — /wings)-^Hi represents the Hi mass contained in the central 
peak of the superprofile. The /cold variable represents the fraction of Hi gas that is in the cold 
phase, which has a narrower velocity dispersion compared to the average turbulent component. The 
velocity dispersion of this gas is likely thermal, and therefore cold H i should not be included in our 



calculation of turbulent energy. Typical measured fractions are /cold ^ 10 — 20% (e.g.. Young et al. 



2003 Warren et al. 2012|). Without strong constraints on individual values for cold Hi fractions. 



we choose /cold = 0.15 for our analysis. 
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6.3.2. Energy from Gravitational Instabilities 



The energy provided by gravitational instabilities over a timescale r is: 



E, 



grav 



Egrav (1 — /wings) (1 — /cold)AfHI T ■ 



(17) 



As before, the (1 — /wings) (1 — /coid)-^Hi factor represents the approximate Hi mass contained in 
the turbulent central peak. 

We now calculate the amount of energy released into the ISM based on the gravitational 



instability proposed by Wada et al. (2002). This instability allows the gas to extract energy from 



rotation instead of from shear, so it is a potential source of energy for dwarf galaxies that lie 



primarily in the regime of solid body rotation. Wada et al. (2002) approximate the energy supply 
rate per unit mass as: 



'grav 



km^ ! 



7^ 



Sgas 



lOMopc^ J VlOOpc 



A 



100 pc 



(18) 



Here, is the angular velocity, Sgas is the gas surface density, A is the scale length of turbulence, 
and hz is the scale height of the disk. We approximate Vt ~ (w2o/2)/(1.5r25), where W2q has been 
corrected for inclination and 1.5r25 is the approximate extent of Hi in galaxies. This approximation 



provides an order-of-magnitude estimate for the angular velocity of gas in our sample (Giovanelli 



k, Haynes 1988 Swaters et al. 2009). We use the (Sgas) measurement for the gas surface density 
and choose A = 100 pc following Tamburro et al. (2009). Based on studies of the scale height in 



dwarf galaxies by Banerjee et al. (2011), we approximate the scale height of our sample as 500 pc. 



We must then find the timescale over which the H i can dissipate its energy. If the H i in the 



central peak is turbulent, this timescale is the turbulent timescale as given by Mac Low (1999): 



TD ~ 9.8 Myr 



A 



lOOpcy VlOkms 



(19) 



where A is the turbulent driving scale and a is the Hi velocity dispersion. Following Tamburro 
et al. (2009), we estimate A ~ 100 pc. We also set a = fXcentrai as measured from the superprofile 
for each galaxy. For our measured cJcentrai values, the turbulent timescale ranges from 9 — 16 Myr. 

To convert the available energy to Hi kinetic energy, the conversion efficiency, egrav = -E'niZ-E'grav 
must be taken into account. We do not adopt a single value for egrav, and instead measure the range 
of the range of egrav that is compatible with our data. In general, implied efficiencies of egrav > 1 
are unphysical in that the Hi component has more turbulent kinetic energy than can be provided 
over the timescale. 



In Figure 19, we compare i^grav to the kinetic energy contained in the central peak. In the left 
panel, we plot -Egrav versus the energy in the central peak. We have shown the background in grey 
to represent the fact that the required efficiencies are > 1 and are therefore unphysical. Dashed 
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lines represent constant Cgrav The left panel shows the inferred value of egrav necessary to drive 
turbulence. 

In all cases, the gravitational instability cannot provide enough energy to drive the observed 
levels of Hi turbulence, falling short by a factor of 10 — 10'^ in spite of the fact that the correlation 
between (Shi) and cxcentrai may have initially pointed at this driver. The discrepancy between 
-E'HLcentrai and -Egrav appears to be more extreme in small galaxies with low values of -Egrav We note 
that in nearly all cases, changing any of our assumptions by a factor of two does not alter the result 
that the Hi harbors far more energy on average across the disk than the gravitational instability 
can provide. The inability of this instability to drive turbulence has been noted before (e.g., 
MacLow|[2004l ), so our results confirm this idea. It is possible that other gravitational instabilities 



are operating to produce the observed (Tcentrai values, but any candidate instability must be able to 
function efficiently in galaxies with low internal shear. 



6.3.3. Energy from Star Formation 



We next turn to star formation as a driver of Hi velocity dispersion. Even though there is 
no straightforward correlation between measures of star formation and (Jcentrai) this assessment 
provides a limit on the efficiencies necessary to couple energy from star formation to the H i gas if 
it is indeed, as widely regarded, the driver of H i turbulence for the bulk of the gas. 

From the measured star formation rate, we can estimate the energy released into the ISM by 
SNe over the turbulent dissipation timescale td as: 



EsF ~ EsF T D 



(20) 



where Esftd is the total amount of energy released from the stellar population over one turbulent 
timescale. This equation implicitly assumes that the rate of energy input from the stellar population 
has been constant over to- 



Tamburro et al. 



(2009), assuming that the 



To estimate Esf , we use the formalism proposed by 
majority of the star formation energy is released by SNe explosions. On average, each explosion 
provides 10^^ ergs of mechanical energy. We take the number of SN per unit stellar mass formed 
to be r/sN = 1-3 x lO^^ SN Mq 



assuming a 



Kroupa (2001 ) IMF with an upper mass limit of 120 



Mq. We next assume that the SFR measured by FUV + 24/im observations has been constant 
over the r^) timescale. The total energy available to the H i due to SNe is then: 



Esf = r/sN (SFR x td) W^'ergs 



(21) 



Since the measured SFRs of our sample are averaged over 10 — 100 Myr, this equation assumes that 
this average measurement is representative of the SFR over td for our galaxies, which is ~ 10 Myr 



(e.g.. Equation 19). This assumption may fail in the case of galaxies with bursty star formation 
histories, since a large recent burst may be able to affect the Hi gas while an older burst may 
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not. However, the energy from supernovae is smoothed out over ~ 50 Myr after a burst, so minor 
fluctuations in the star formation history may not be too important. It also assumes that the initial 
mass function (IMF) is well-sampled, which may be an issue for galaxies with very low SFRs. This 
formalism indicates that the galaxies with the lowest SFRs (~ 0.5 — 1 x 10~^ yr~^) would 
have approximate 1 SNe over 1 — 2 x 10^ yr. We note that stellar winds can provide additional 
mechanical energy into the ISM. 

As with -Egrav) the conversion efficiency between energy available from star formation and Hi 
kinetic energy, esF = Em/EsF, must be taken into account. Values of esF > 1 are unphysical, 
as in these cases star formation provides less energy than is contained in the central Hi peak. 



Additional limits have been placed on esp by simulations. Thornton et al. (1998) found that the 
average efficiency (esp) ~ 0.1, while other simulations measure efficiencies that can be as high as 



0.5 ( |Tenorio-Tagle et al.|[T99l| ). As with 

£grav; we measure the range of esF that is compatible with 



our data. 



In Figure 20, we compare the energy provided by SNe over one turbulent timescale to the 
energy in the Hi gas. The format is the same as Figure [T9j In the left panel, we plot the star 
formation energy versus the kinetic energy in the central Hi peak. The dashed grey lines show 
constant esp- The unphysical region where esp > 1 is shown in dark grey. The simulations by 



Thornton et al. (1998) suggest that esp is never higher than 0.1; we shade regions above this 



threshold in light grey. The left panel shows the inferred values of egp necessary to drive turbulence 
with star formation energy. Again, we have shown efficiencies above the more stringent theoretical 



maximum of 0.1 from Thornton et al. (1998) in light grey. 



For all galaxies, we find that recent star formation provides enough energy to drive the observed 
turbulence over a single turbulent timescale with efficiencies 0.01 < esp ^ 0.1, well in line with 



the limits from both Tenorio-Tagle et al. (1991) and Thornton et al. (1998). However, it is likely 



that some fraction of this energy (~ 35%, e.g., Joung et al. 2009) goes into accelerating Hi from 
its undisturbed state into the wings of the superprofile, as explored in the following section. In 
this sense, these efficiency estimates are lower limits. On the other hand, some of the mass in the 
central peak is likely kinematically associated with the wings, thus lowering esp values. 



6.3.4- Superprofile Wings 

Next we compare the kinetic energy contained in the wings of the superprofiles to that provided 
by star formation. Although the wings contain only a small fraction of the H i mass ((/wings) ~ 0.11), 
the velocities are very high (((Twings/ccentrai) ~ 3), suggesting that they may harbor a significant 
fraction of the kinetic energy. 

We calculate the kinetic energy in the wings as follows. If all Hi in the galaxy started in the 
central turbulent component, some extra energy is necessary to accelerate the gas from v ~ fJcentrai 
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to V 



wmgs • 



The excess kinetic energy in the wings is then: 



E, 



wmgs 



fmngs -^HI ((^wings '^central) • 



(22) 



Based on our definition of ciwings (Equation 10), this gives the total energy necessary to accelerate 



a mass /wingsAfni from a Gaussian velocity distribution into the observed wings of the superprofile. 
Because some of the Hi in the central peak may be kinematically associated with the wings, this 
assumption provides a lower limit on the energy contained in the wings. 

We next must choose a timescale over which to consider energy input from star formation for 
the wings. The relevant timescale to consider is not necessarily straightforward because the source 
of kinematics in the wings is unclear. If the wings are representative of bulk motions, such as away 
from star forming regions, the relevant timescale should be related to how long these bulk motions 
are expected to persist. On the other hand, if the wings are turbulent, the best timescale may be 



the turbulent timescale given in Equation 19 Because the source of the wings is not necessarily 
clear and because neither case provides a definitive timescale, we assess the ability of star formation 
to provide enough energy over both timescales in turn. 

First, we consider the scenario where the wings represent bulk gas motion away from star 
forming regions. To estimate the relevant timescale associated with this component, we capitalize 
on recent studies of Hi holes, as these structures often exhibit velocity structures similar to what 
we expect to find in the wings. The kinematic age of H i holes can be estimated based on their size 
and expansion velocities, if observable. This calculation is very uncertain and typically provides 
at most an upper limit to the true age as expansion velocities are expected to slow over time. 
However, they provide an order-of-magnitude estimate of the timescale over which these structures 



are observable. Seven galaxies overlap between our sample and that of Bagetakos et al. (2011) 



(NGC 2366, Holmberg II, IC 2574, Holmberg I, NGC 4214, DDO 154, and NGC 7793), who find 
a mean kinematic age in dwarfs of ~ 32.5 Myr. We therefore adopt this value for Thoies with the 
caveat that it is uncertain, and only consider star formation energy input over the approximate 
timescale on which H i signatures of expansions are expected to decay. We use the same formalism 



as in Equation 21 and § 6.3.3 to estimate star formation energy, substituting 32.5 Myr for r/j. 



In Figure 21 , we compare the energy available from star formation to the energy in the super- 
profile wings. In the left panel, we plot the star formation energy versus the wing kinetic energy for 
each galaxy. Dashed grey lines indicate a constant efficiency of transferring star formation energy 
to kinetic Hi energy, where esF,wings = -E'wings/-£'SF- As in Figures 19 and 20, grey regions of the 



plot show unphysical or theoretically prohibited efficiencies, i.e., esF.wings > 1 and esF,wings > 0.1. 
The right panel shows the inferred efficiency esF, wings versus the star formation energy available. 

We find that all of galaxies require efficiencies of only < 0.05 to produce enough energy to drive 
H I gas into the wings over ci timesccilc of 32.5 Myr. The distribution licis a niGdictn tgp^wings 

= 0.013 

with a standard deviation of 0.007. These estimates are well below with the theoretical maximum 
of 0.1 — 0.5 found by Tenorio-Tagle et al. (1991) and Thornton et al. (1998), and are in line with 
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simulations by Joung et al. (2009) that show that ~ 35% of star formation energy goes into driving 



large-scale bulk motions. The required efficiencies are also much lower than estimates of efficiencies 



necessary to drive larger Hi holes, which range between 1 - 40% at their kinematic ages (e.g., Weisz 
et al. 2009; Warren et al. 2011; Cannon et al. 2011 Bagetakos et al. 2011). However, we have 



derived these estimates from global properties. A more precise determination is necessary using 
spatially-resolved data scales, as previous studies have shown that the value of the second moment 



declines with radius Tamburro et al. (2009), therefore changing the energy in the wings based on 
location in the galaxy. 

We next consider the possibility that the gas in the wings is instead representative of a turbulent 
component. In this case, the relevant timescale to consider is the turbulent dissipation timescale, 



as given by Equation 19 Even though the velocity profile of H I in the wings is not Gaussian in our 
parameterization, it again provides an order-of-magnitude estimate of how long this component can 
dissipate its energy if it is indeed turbulent. In this case, it may be more relevant to use measured 
CTfc values for double Gaussian fits instead of cJwings, because our parameterization explicitly removes 
gas with small velocities from the Cwings calculation. Because the double Gaussian fits appear to 
be determined primarily by the wings (see Appendix O), we choose to use erf, to calculate td as 



the second timescale to consider. If we substitute Uh for a in Equation [T9| we obtain values for 
r£) ~ 6 — 10 Myr, or ~ 70% of those determined for the central component. As before, we use 



Equation [21] to calculate the input energy over this turbulent dissipation timescale determined for 
the wings. 



Figure 22 shows the comparison of energy in the wings to energy provided by star formation 



over one turbulent timescale. The format is the same as Figure 20 In this case, implied efficiencies 
are higher by a factor of ~ 3 — 5 due to the difference in timescales, and are much closer to the 



theoretical maximum of 0.1 found by Thornton et al. (1998). 



6.3.5. The Entire Superprofile 

In this section we assess whether star formation provides enough energy to produce the H I ve- 
locity distribution seen in the superprofiles. In this case, the kinetic energy of the entire superprofile 
is simply: 

EsP = ^^M{v)v^ (23) 

V 

where M{v) is the Hi mass at velocity v. 

As with the central peak and the wings, we must determine the timescale over which to 
calculate energy input from star formation. Since the relevant timescales are unclear, we use both 



the turbulent timescale for the central peak and a fixed timescale of r = 32.5 Myr in Equation 21 
to calculate energy input. 



We show the comparison between energies over the turbulent timescale in Figure 



23 



Over this 



- 41 - 



timescale, star formation provides enough energy to drive the fuU shape of the superprofile, but 
many galaxies he in the e > 0.1 region unfavored by simulations (Thornton et al. 1998). 



Figure 24 again shows the comparison between star formation and Hi energies, but over the 
the fixed timescale of r = 32.5 Myr. Over this timescale, star formation provides enough energy 
to drive kinematics in the entire Hi superprofile at with e < 0.1. We note, however, that it is 
likely that the central peak and the wings of the superprofile have different associated timescales 
and efficiencies, so choosing single values to represent both components may not be physically 
appropriate. 



6.4. Estimating the Scale Height of Hi in the Sample 



The scale height of the Hi layer perpendicular to an isothermal, self-gravitating disk can be 
determined based on its velocity dispersion and disk surface mass density ( van der Kruit|1981 ). We 



use a method similar to that presented in Ott et al. (2001 ) and Warren et al. (2011 ) to approximate 



Hi scale heights for our sample. The scale height hz is given by van der Kruit (1981) as: 



'gas 



(24) 



where (Xgas is the velocity dispersion perpendicular to the disk, G is the gravitational constant, and 
Pt is the stellar mass density of the disk. If we assume that H i has a Gaussian distribution in the 
z-direction, we find that: 

A'^Hi = \/27r/i^nHi,o, (25) 
where A'^hi is the H i column density and nni o is the number density at the midplane of the disk. 



These two equations were combined by Ott et al. (2001) to find an expression for Hi scale height 
in terms of observables: 



hz = 579 



gas 



lOkms 



10^^ cm ^ 



PHI 

Pt 



pc. 



(26) 



where pm/pt is the ratio between Hi density to total disk density, and can be approximated as 
{pm/pt) = (Mhi/Mt) where Mt = 1.36Mhi + M^. 

Using our superprofiles and measured galaxy properties, we can estimate the average H i scale 
height of the disk. We first assume that the velocity dispersion is isotropic, such that Cz — (^central* 
Second, we convert (Shi) to A'^hi units with the caveat that we have averaged these quantities 
over the disk and have not included the contribution of dark matter. We note that since these 
values are averaged over the disk of the galaxy, any spatial variation in these parameters is no 
longer distinct. Hi scale heights in galaxies are expected to fiare at large radii, as gas velocity 
dispersions remain relatively constant but disk surface density declines with radius. This method 
therefore gives an estimate for H i scale height that is weighted toward the regions with the highest 
Hi surface densities. The average Hi scale heights derived from these values are listed in Table |8| 
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The majority of galaxies have imphed scale heights between 100 < {hz 
of 320 pc and interquartile range of 210 — 480 pc. 



< 700 pc, with a median 



Scale heights for some of the galaxies in our sample have been determined using other methods. 



Banerjee et al. (2011) modeled the dark matter and baryonic components of halos for DDO 154, 
Ho II, IC 2574, and NGC 2366. They obtained scale heights between ~ 130 pc at r = kpc to 
~ 1 kpc at r = 6 kpc for DDO 154; between ~ 180 pc at r = 1 kpc to ~ 1 kpc at r = 7 kpc for 
NGC 2366; between ~ 350 pc at r = 1.5 kpc to ~ 700 pc at r = 9 kpc for IC 2574; and a fixed 



scale height of ~ 400 pc at all radii for Ho II. Compared to the Banerjee et al. (2011 ) scale heights. 



our method yields scale heights that are within the same range for DDO 154 and NGC 2366 but 
smaller by 35% and 25% for IC 2574 and Ho II, respectively. 



6.5. A Universal Hi Profile Shape? 



As seen in Figure [Sj the superprofiles show a distinct velocity distribution: a central peak 
with strong contributions from non-Gaussian wings. After normalization to the same HWHM, the 
uniformity of the profiles is striking, especially considering that the low-mass dwarf galaxies in our 
sample are typically characterized by irregular velocity fields and morphologies, and stochastic. 



varied star formation histories (Weisz et al. 2011). The residuals also show a surprisingly similar 
shape, with the peak often occurring at 2 x HWHM across the observed range of /wings- This global 
regularity exists in spite of the fact that the individual line-of-sight H i profiles have a much more 
varied shape, with some showing asymmetry (e.g.. Young et al.||2003 Warren et al.||20T2 ) or double 



peaks indicative of expanding structures (e.g., [Bagetakos et al. 2011). Statistically, however, the 



sum of these profiles generates the same kinematic distribution from galaxy to galaxy. 

The shape of the superprofiles on global scales is qualitatively similar to those found by other 
studies of average Hi line profiles, which also show a mostly Gaussian central peak with broader 



wings (e.g., Dickey et al. 1990; Boulanger & Viallefond 1992 Kamphuis & Sancisi 1993 Braun 



1997 Petric &: Rupen||2007 lanjamasimanana et al.|2012 ). It also matches the shape of simulated 



Hi profiles found by Joung et al. (2009) for a supernova-driven turbulent medium, though at lower 



star formation intensities. However, quantitative comparisons among the surveys are hindered by 
the vast differences both in observational parameters such as spatial and spectral resolution as well 
as in techniques for removing the rotational velocity. Nonetheless, many of the studies, including 
our own, find that the broad central peaks can be characterized reasonably well by Gaussians with 
widths of 5-10 km s~^, with additional wings to either side. The similarities indicate that the shape 
of the H I line profile, and therefore the general kinematic structure of H i line profiles, are relatively 
independent of galaxy properties. However, uniform studies of larger spirals with better velocity 
resolution, now possible with the larger bandwidth of the newly-updated JVLA, are necessary to 
confirm this idea. 



7. Conclusions 



We have generated a measure of global Hi kinematics in a sample of nearby dwarf galaxies 
from VLA-ANGST and THINGS by summing the contribution to a global line profile for each 
line-of-sight spectrum after removing rotation from each spectrum. The resulting superprofile for 
an individual galaxy provides an intensity- weighted average of its individual H i line profiles. 

We interpret the superprofiles as composed of a central peak indicating average turbulence 
with higher-velocity wings to either side. We parameterized them with four parameters describing 
the width of the central peak (cJcentrai), the characteristic velocity of the wings ((T„ings), the fraction 
of gas in the wings (/wings), and the asymmetry (a). We have compared these parameters to 
various global galaxy properties in order to determine what, if any, physical causes are behind Hi 
kinematics. 



• The dynamic range of Cccntral is quite small, varying only between ~ 6 — 10 km s~ across our 
sample. We find a correlation between fXcentrai and (Shi) in the dwarf sample which is not 
significant once higher mass galaxies are added. The measured cicentrai values are close to but 
slightly higher than line widths that can be driven by background UV heating. It is possible 
that base H i velocity dispersions are set by this heating, with star formation imparting only 
additional energy. 

• The characteristic velocity of gas in the wings, (Twings, increases with (Ssfr), (Sbaryon); and 
(Shi), implying that star formation could be one way to accelerate Hi to velocities faster 
than expected compared to the surrounding turbulent medium. 

• The fraction of gas in the wings, /wings, increases with galaxy mass and with SFR / Mhi, so 
galaxies with relatively high SFR or low Mhi could be better able to accelerate H i to higher 
velocities. 

• The asymmetry, a, decreases with both SFR and with galaxy mass, and is primarily in the 
wing regions. This supports the idea that star formation can accelerate Hi away asymmetri- 
cally, so galaxies with smaller SFR likely have fewer star-forming regions and thus show more 
asymmetry. 



In all cases, our trends exhibit large scatter. Since many of the physical properties we exam- 
ined vary on both radial and spatial scales, future analyses must incorporate this information to 
disentangle any causal connection with Hi gas kinematics. 

We have also compared the energy contained in the H i superprofiles with the energy provided 



by the gravitational instability from Wada et al. (2002) and by star formation. We find that this 



gravitational instability cannot provide enough energy to drive turbulent line widths on timescales 
of ~ 10 Myr, while star formation can. Star formation also imparts enough energy to accelerate 
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gas into the wings of the profile over timescales of ~ 32.5 Myr, with imphed efficiencies below the 
theoretical maximum of 0.1 - 0.5. 

We derived average H i scale heights for the sample, with most galaxies exhibiting scale heights 
of a few hundred pc. 

Finally, we found that the average Hi superprofile shape, when scaled to the same HWHM, 
has a remarkably similar shape from galaxy to galaxy, with variations primarily in the wings at 
low levels. The shape of the central component differs from a Gaussian by only ~ 0.05 with more 
varied wings showing additional emission 5-10% above the scaled Gaussian fit. This similarity 
implies that the physical processes setting the kinematics of H i in galaxies function similarly in all 
dwarf galaxies. 
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Science foundation operated under cooperative agreement by Associated Universities, Inc. Support 
for this work was provided by the National Science Foundation collaborative research grant "Star 
Formation, Feedback, and the ISM: Time Resolved Constraints from a Large VLA Survey of Nearby 
Galaxies," grant number AST-0807710. This material is based on work supported by the National 
Science Foundation under grant No. DGE-0718124 as awarded to A. M.S. 



A. Noise Estimates for the Superprofiles 

The fiux rescaling process is important for interferometric data because it provides an accurate 
measurement of the true Hi fiux of each galaxy. Since this process only rescales the residuals, the 
noise properties of the rescaled data cube are complicated. While the highest noise peaks can 
show up in the clean components of the rescaled data cube, the majority of the noise is still in 
the residuals. Traditional estimates of noise therefore provide an overestimate of the noise for 
superprofiles generated from rescaled data cubes. 

The noise on a single point of a superprofile generated from a standard cube can be approxi- 
mated as: 

0"SP = CTchan X ^J ^^pix/A'pix/beann (Al) 

where fXchan is the rms noise level in a single channel of the data cube, A^pix is the number of pixels 



contributing to each superprofile point, and A^pix/beam is the number of pixels per beam (12012). 



In the left panel of Figure 25 we apply this formula to the superprofile generated from the 
standard cube and to that from the rescaled cube. The upper panel shows the absolute fiux 
measured in each superprofile. In the lower panel we have normalized the superprofiles so that the 
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maxima are the same. It is clear that the fractional noise is much larger in the rescaled cube. 

In the right panel, we have rescaled the noise for the superprofile generated from the rescaled 
cube by the ratio of total fluxes, i^rescaied /-^standard such that the noise is given by Equation |8j 
This rescaling produces a noise estimate with a similar fractional uncertainty compared to the 
standard cube. This estimate is not exact, because the highest noise spikes are contained in the 
clean components. However, the majority of the noise is in the residuals, which are rescaled. 
This method also provides a noise estimate that matches the fractional uncertainty on each point 
compared to superprofiles generated from the standard cubes. Since it is a better representation 
of the fractional uncertainty on each superprofile point, we adopt Equation |8] when calculating the 
effects of noise on the superprofiles. 



B. Effects of Observational Settings on Superprofile Shapes and Measured 

Parameters 



We have performed a number of tests to ensure the validity of our results and to estimate 
uncertainties on the measured parameters, cJccntrah o'wings, /wingsi and a. In particular, we have 



examined the effects of fpeak uncertainties (B.l); finite spatial resolution (B.2); finite velocity res- 



olution (B.3); and noise on each superprofile point (B.4). Finally, we review the final uncertainties 



on the measured parameters (B.5). 



B.l. Uncertainties in Vp^ak 

We explore how our superprofile parameters are affected by the uncertainties in determining 
Vpeak, which could possibly generate broader superprofiles or more fiux in the wings. 

For each of our four test galaxies (GR 8, Sextans A, UGC 4483, and NGC 2366), we start with 
the pixel-by-pixel uncertainty in determining "^peak as a function of S/N ratio, determined from our 
Monte Carlo tests for the four test galaxies (§ 3.2, Figure [6]). We first assume that all pixels would 



contribute a Gaussian with a width of (Jcentrai in the absence of any uncertainties on Upeak- For each 
pixel above our S/N > 5 threshold, we generate a random offset drawn from a Gaussian with that 
pixel's standard deviation of determining t^peak and with the S/N as the amplitude. We generate a 
fake superprofile by summing all of the (Tccntrai Gaussians with their respective velocity offsets. 



The results of this test are shown in Figure 26 for our four test galaxies. For each of the four 



galaxies, the upper panel shows the "observed" fake superprofile as the black solid line; the input 
superprofile we would have expected in the absence of any uncertainties on fpeak as the blue dash- 
dot line; and the HWHM-scaled Gaussian fit as the dashed red line. In all cases, the differences 
are smaller than the line widths. The lower panel shows the residuals (i.e., "observed" - fit and 
"observed" - input). 
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The differences between the input Gaussian, tfie "observed" fake superprofile, and the fit are 
< 0.005 in all cases. We also find that the width of the superprofile is increased by < 0.5%. The 
uncertainties in I'peak therefore have a negligible effect on the superprofile shapes and parameters. 



B.2. Finite Spatial Resolution 



The combination of finite spatial resolution and rising rotation curves at the centers of galaxies 
can increase the width of observed H i line-of-sight spectra in the central regions, which could then 
either increase the width or mimic Hi in the wings of the observed superprofile. Our sample 
of dwarf galaxies likely have either slowly-rising rotation curves (DDO 154 and NGC 2366; e.g., 



de Blok et al. 2008) or primarily display solid body rotation typical of dwarfs (Oh et al. 2011). 



Because these rotation curves have a smaller gradient with increasing radius, we expect to see less 
of an effect from beam smearing in the central regions compared to larger spiral galaxies, but we 
must still to understand its effects. 

To quantify the effects of beam smearing, we have developed a Python module to generate a 
suite of model galaxy observations using NGC 2366 as our test galaxy. This galaxy has the steepest 
rotation curve of our sample, and thus would be the most affected by this particular bias. We use 
the observed Hi surface brightness distribution plus the observed inclination and position angle 



from de Blok et al. (2008) to generate the model Hi distribution. We also assume that the disk has 
an exponential distribution in the z-direction with a scale height = 500 pc, a typical observed 



value for Hi scale heights in dwarfs (e.g., Banerjee et al. 2011 Warren et al. 2011). Changing this 



assumption to either 100 pc or 1 kpc does not strongly influence our results. 

We next impose a rotation curve that can be modeled as a linear rise for radii smaller than 
Tflat, with a flat regime at larger radii with circular velocity vug^t- 

Finally, we assume that all line-of-sight H i spectra have Gaussian velocity distributions with 
a dispersion of 6 km s~^. While this assumption is not necessarily indicative of the true dispersion 
as a function of radius (e.g., Tamburro et al.||2009 >, it does allow us to quantify the effects of beam 
smearing on a uniform H I profile. In order to estimate the effects of declining velocity dispersion, 
we also generate a model where the intrinsic velocity dispersion is chosen by an exponential fit to 
the radial average of the second moment map. 



For this test, we use three models: one with the observed rgat = 1-9 kpc, i^fjat = 60 km s~^ 
(de Blok et al. 2008), and cthi = 6 km s~^; a second with an extreme r^at = 0.5, t^flat = 60, and, 
Chi = 6 km s~^ ; and a third with rgat = 1.9 Vfiat 
the second moment map. 



60, and uhi (r) set by the exponential fit to 



To place H i clouds in the model cube, we draw a sample of points from the H i surface brightness 
distribution, assuming each cloud represents a gas cloud at the observed spatial position. We then 
distribute these points randomly in the z-direction using our assumed exponential z distribution. 
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Each point is smoothed in velocity space into a Gaussian with the central velocity determined by 
the rotation curve at that radius and the width dependent on the model. Finally, we scale the cube 
so that the total Hi mass is the same as measured in the galaxy. This yields a cube that has Hi 
line-of-sight spectra unaffected by spatial resolution with a velocity resolution of 2.6 km s~^. 

To reproduce the effects of finite spatial resolution, we smooth the cube to our working 200 
pc resolution using a circular Gaussian beam with a FWHM of 200 pc. For both the true and 
convolved cube, we find Vpeak for each pixel and then generate a superprofit using the same method 
as described in f|3j We fit and parameterize the superprofit using a single Gaussian scaled by the 
superprofit's HWHM (Q). 



We show the resulting true and smoothed model superprofits compared to the observed su- 
perprofit in Figure [27} In all cases, the superprofit from the smoothed cube shows very small 
differences compared to the true cube. The measured Gaussian dispersion is slightly wider by 
Aucentrai, spatial ^ 0.5 km s~^. Additionally, the convolved cubes have a negligible fraction in the 
wings (A/„ings_spatiai ^ 0.01), a value much lower than the typical range observed in our sample 
(0.05 < /wings < 0.15). With such a small contribution to the wings, Swings should not be strongly 
affected by beam smearing. We note that the superprofit generated from the model with exponen- 



tially declining am is well-fit by a Gaussian with a width ~ 10 km s ^, similar to the Tamburro 



et al. (2009) results. None of these models is able to reproduce wings at the observed magnitude. 



In summary, while finite spatial resolution does contribute a small amount of broadening, it is 
at a low level compared to the observed widths and is not strong enough to generate spurious fiux 
in the wings. 



B.3. Finite Velocity Resolution 

We also examine the effects that finite velocity resolution has on our results. Since our sample 
is composed of galaxies observed with a variety of velocity resolutions (0.6, 1.3, and 2.6 km s~^), 
we must quantify any effects that arise from these differences. For each galaxy, we bin the observed 
standard and fiux-rescatd cubes to the coarser velocity resolutions of our sample. We then find 
Vpeak, generate a superprofit, and measure parameters for each new resolution. 

In Figure [28] we show the superprofits generated from the binned cubes for DDO 125. The 
upper panels shows the superprofits themselves, and the lower panels show the difference between 
superprofits generated from the binned cube and the original cube. The superprofits behave 
similarly for Av = 0.6,1.3 and 2.6 km s~^, but the superprofit for the Av = 5.2 km cube 
shows > 5 % differences compared to the original. Similar behavior is evident in all the superprofits 
generated from cubes binned to Av = 5.2 km s~^, and thus galaxies whose original data cubes have 
Av = 5.2 km s~^ have been excluded from our sample. 



We show the resulting parameters as a function of bin size for all galaxies in Figure 



29 



The 
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four panels in the plot are cTcentrah Swings; /wings, and a. Each color represents a different galaxy. We 
find that three of our four parameters (ccentrai; ""wings 

and /wings) are relatively well-behaved with 
increasing velocity resolution, showing minor variations relative to the range seen in the sample. 
However, the standard deviation of variations in a is ~ 0.06; we account for this variation in the 
final uncertainties ( ^B.5 ). 

In general, coarser velocity resolution slightly increases the width of the central component 
and the velocity of the wings, which then places slightly less flux in the wings. Comparing the 
same galaxy at 1.3 km s^^ to 0.6 km s~^, we find median differences of Acrcentrai,vei = 0.08, 
Ao"„ings,vci = 0.13, and A/^ings.vei = —0.005. Comparing between 2.6 and 1.3 km s~^, we find 
median differences of A(Tccntrai,vci = 0.17, Acj„ings,vei = 0.32, and A/^ings.vei = —0.006. 



B.4. Noise 

The final influence on the measured parameters is the noise of the superprofits and is especially 
important for the lowest-mass dwarfs. To quantify the effects of noise on the measured parameters. 



we start with the noise estimate for each point (Equation [8| § 3.4). For each sample galaxy, we 
assume that the measured superprofile is true. We then add noise to each point drawn from a 
Gaussian distribution with a width usp calculated using Equation [8} Finally, we measure the 
parameters from this "noisy" superprofile. 

After repeating the above process 10,000 times, we have obtained an estimate of the typical 
range of parameters allowed by the noise on the superprofiles. For each parameter, we fit a Gaus- 
sian to the histogram of "noisy" parameter values and adopt its width as the uncertainty due to 
superprofile noise on each parameter. 

The median uncertainties and interquartile ranges on the parameters due to noise are: Adcentrai.n 

^■'^^oioi S ^"^wings, noise — f-9o!8 ^ ^/wings, noise — O-'-'l^oioOEi ' and Aanoisc — 0-0"^Oo'o09 



B.5. Overview of Parameter Uncertainties 

We use the above tests to estimate the total uncertainty of the measured superprofile parame- 
ters. We include only uncertainties that have a non-negligible effect on each parameter in its total 
uncertainty estimate by adding the uncertainties due to various observational effects in quadrature. 

For each galaxy, the final uncertainty for ^central is^ 

Adcentral = ^ (0.5 km + (0.17 km + (Aacentral,noise)2. (Bl) 

The 0.5 km s~^ and 0.17 km s~^ errors are due to the effects of finite spatial resolution and finite 
velocity resolution. The value for Adcentrai, noise is different for each galaxy. We do not include errors 
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from Vpeak uncertainties, as they are an order of magnitude lower than the other uncertainties. The 
^o'ccntrai valucs are dominated by errors due to finite spatial resolution. 

The final uncertainty of fiwings is: 

AcTwings = \/ (0.13 km S"^)2 + (A(J„ings,noisc)^- (B2) 

The 0.13 km s~^ error is due to the effects of finite velocity resolution, and Acrwings,noise is different 
for every galaxy. We do not include errors from fpeak uncertainties or finite spatial resolution, as 
the small amounts of flux they contribute to the wings do not strongly clicinge cr^jj^gg values. The 
Auwings values are dominated by uncertainties due to noise on the superprofiles. 

The final uncertainties for /wings is: 

A/wings = ^ (0.01)2 + (AU-^g 

s, noise 

)2. (B3) 

The 0.01 uncertainty is due to finite spatial resolution. As before, A/„ings,noise is different for each 
galaxy. We neglect uncertainties due to finite velocity resolution and Wpeak uncertainties, as they 
are over an order of magnitude smaller than the included uncertainties. Both values contribute 
approximately equally to the final uncertainty A/„ings. 

The final uncertainty for a is: 

Aa = V(0.06)2 + (AOnoise)2. (B4) 

The 0.06 uncertainty is due to finite velocity resolution, as the peak determination is less precise 
in cubes that have larger velocity resolutions. We do not include uncertainties from finite spatial 
resolution in this estimate. The included uncertainties contribute fractionally different amounts in 
different galaxies, so Aa is not typically dominated by either uncertainty. 



C. The Interpretation of Double Gaussian Fits 

In a recent study, |I2012| parameterized superprofiles for THINGS galaxies using double Gaus- 
sian fits and then argued that the two components were representative of cold and warm H i in the 
galaxies. We had independently pursued this approach for our study based on the low reduced 
value of double Gaussian fits compared to single Gaussian fits. Unlike [120 12 we chose to abandon 



this parameterization in favor of the simpler HWHM parameterization for reasons described in this 
appendix. 



We have previously given results for our double Gaussian fits in § 4.3 and Table [5] However 



a number of potential concerns about the physical meaning of the double Gaussian fits arose as we 
explored them in more detail. 

We first refer to the strong similarity in shape seen in Figure |8| when all superprofiles are scaled 
by their HWHM value. We now question whether the same similarity is seen when scaling by width 
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of the narrow or broad Gaussian components, which might be expected if (T„ or ab were physically 
relevant quantities. In Figure [30] we again show the scaled superprofiles for all the galaxies, but 
we now scale the velocity axis by a different measured parameter in each panel: HWHM (upper 
panel); the width of the narrow component of the double Gaussian fit, an (middle panel); and 
the width of the broad component of the double Gaussian fit, ab (lower panel). The left column 
shows the superprofiles after normalizing the scaled velocity axes so that the median width of the 
superprofiles is aligned, to better show the variation in shape. We then plot the median superprofile 
in red. The right column of Figure |30] shows each superprofile minus the median superprofile for 
that scaling over the same velocity range. 

The HWHM scaling removes most of the variations in shape among the profiles. In contrast, the 
specific values of cj„, and ab have little direct bearing on the overall profile shape. Quantitatively, the 
rms residuals around the median scaled superprofile are 5 times larger for the superprofiles scaled 
by either an or at compared to those for the HWHM-scaled superprofile. The HWHM-scaling 
appears to provide the best characterization of the superprofile shape. 

We next ask if the Gaussian components behave similarly when scaled to the superprofile 
HWHM. Since the HWHM-scaled superprofile shapes are very similar, we would expect the narrow 
and the broad components of the double Gaussian fits to have similar properties when scaled in 



the same way. We plot the results of this test in Figure 31 The left panel shows the superprofile 
for each galaxy scaled to its HWHM. We then overplot the two Gaussian components, which have 
also been scaled to the same HWHM as the corresponding superprofile. The upper panel shows the 
narrow component, while the lower panel shows the broad component. At first glance, the shapes 
of the scaled Gaussian components are strikingly different compared to the similarity in the global 
superprofiles' shapes. In other words, the properties of the two Gaussian components vary wildly, 
while somehow conspiring to preserve the same overall shape. 



The right panel of Figure 31 is the same as the left panel, but for clarity only Sextans B 
and NGC 7793 are shown. Although the upper ~ 50% of the superprofiles are nearly identical, the 
best-fit double Gaussians are very different. Since the superprofile for NGC 7793 has broader wings 
relative to the superprofile HWHM, the broad Gaussian is forced to a larger ab value and a smaller 
amplitude to fit the wings. The narrow component is then left to make up the remainder of the 
central peak. Because the wings are broad, which leads to a low-amplitude broad Gaussian, the 
narrower component is forced to a high amplitude to match the center of the profile. In contrast, 
the superprofile for Sextans B has lower level wings. The broad Gaussian component therefore 
has a higher amplitude, leaving less of the central peak for the narrow component. This behavior 
indicates that the relative amplitude and fractional area of the broad component are strongly driven 
by the shape of the wings. Even though the wings of NGC 7793 have only ~ 9% more flux than 
those of Sextans B, the relative amplitudes of the broad components differ by nearly a factor of two, 
and the fractional areas of the broad components differ by more than a factor of three. The large 
differences in double Gaussian parameters for superprofiles with remarkably similar shapes is the 
first indication that the double Gaussian fits may not be tracing physically-meaningful properties. 
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We see further evidence for this behavior when we quantitatively explore the correlations among 
the double Gaussian parameters. To provide a relative comparison, we have generated superprofiles 
for galaxies in 12012 following their methodology for their well-behaved "clean" subsample. In 
particular, we have used the natural weighted cubes at their instrumental resolution and included 
all pixels with S/N > 3. We do not apply any masking. We can then calculate the HWHM of each 
superprofile. 



Based on Figure [STj we have speculated that the relative strength of the narrow and broad 
components are being set primarily by the relative velocity of the wings compared to the central 



profile (i.e., crj/HWHM). We explore this idea in Figure 32 where we plot the amplitude ratio 
between the narrow and broad Gaussian components versus Ufe/HWHM. The left panel of the plot 
shows the galaxies in our sample, using the double Gaussian fits to the superprofiles discussed in 
this paper (Table [5]). The right panel shows the equivalent data for the clean subsample of 12012 



For this comparison, we calculate the Narrow / Broad amplitude ratio, based on the formula for 
the area under a Gaussian, as follows: 

Narrow / Broad Amplitude = —r- — , (CI) 

where An/Af, is the ratio of the areas under the narrow and broad components, and (T„ and fif, are 
the widths of the narrow and broad Gaussian components. All numbers are taken directly from 



12012 In this case, we have taken at directly from 12012 and calculated the HWHM from our 



12012-like superprofiles. 

It is immediately clear that profiles with broader wings relative to their characteristic HWHM 
width have narrow components with much higher amplitude ratios compared to the broad compo- 



nent. This correlation is a quantitative representation of the behavior shown in Figure 31 the broad 



Gaussian component primarily fits the wings, leaving the narrow Gaussian to fit the remainder of 
the superprofile as best it can. This behavior is exacerbated by the weighting scheme. Be- 

cause this scheme gives pixels with fewer contributing points more weight, the superprofile wings, 
which are produced by only a small fraction of the Hi, are weighted most strongly. Since the 
relative broadness of the wings strongly affects the relative amplitudes of the two Gaussians, it is 
unclear if the two components are tracing the same type of gas in superprofiles with varying levels 
of wing importance. 

We can also look at the ratio of narrow to broad component areas compared to the fraction 
of gas in the wings of the profile (/wings) as measured in § 4.2 In Figure 33 we plot An /A), versus 
/wings- Again, the left panel shows data from our sample, while the right panel shows that from 



12012 For the latter, we have taken An/Ah directly from 12012, and we have used our 12012-like 
superprofiles to measure /wings with the caveat that we have not brought the 12012 sample to a 



common physical resolution. If the broad Gaussian fit were indeed tracing a warm component, 
then on average we would expect galaxies with more flux in the wings to have a higher fraction 



of gas in the broad component, leading to lower A^/A^ values. For the 12012 numbers, we find 
an unexpected positive correlation, where galaxies with more fiux in the wings, as measured with 
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/wings ) have a higher A„/A{, values and therefore a higher fraction of gas in the narrow component. 
We do not see a similar trend in our data, but it is possibly due to the fact that we are measuring 
smaller galaxies with lower overall S/N ratios, higher asymmetries, and a smaller range in /wings- 

These two figures call the physical interpretation of the double Gaussian fits into question. It is 
unclear if the parameters are measuring physical quantities in the galaxies, or if they simply reflect 
the strength of the wings. The correlations between physical properties and both cin/cTft and A„/Afe 
found by |I2012] were attributed to star formation, which could be responsible for driving gas into 
the wings of the profile. We have discussed a similar idea in this paper using a different superprofile 
parameterization. Surprisingly, though, An/A^ and (Jn/(^b show, if anything, only weak trends with 
the direct measure of star formation represented by Hq luminosities. Additionally, since the wings 
of the profile appear to set the relationship between a number of the other measured parameters, 
it not clear if the narrow component is truly tracing cold Hi or if it is just another reflection of the 
superprofile wings. 

We can also turn to constraints provided by previous studies of individual Hi line-of-sight 
spectra. Many previous studies of spatially-resolved H i line profiles in external galaxies have found 
that the line-of-sight profiles are often well- fit by single Gaussians, with only < 20% of profiles 



exhibiting non-Gaussian structures such as broad wings or asymmetries (e.g., Young et al. 2003 
Warren et aljlmll ). The narrow and broad Gaussians fit to these small number of line profiles are 
often interpreted as the emission from CNM and WNM. However, it is unclear if this interpretation 
extends to the statistical measurement of Hi line shapes measured by the superprofiles, since 
information from each individual line profile is no longer distinct. For example, the asymmetric 
line profiles, when added together into a superprofile, can combine to form broader wings than 
what would be measured on a spatially-resolved basis. The spatially-resolved studies of Hi line 
profiles have estimated the fraction of Hi mass in the cold component to be only ~ 20% of the 
total Hi mass, which implies A„/A5 values of ~ 25% - much less than those measured by double 
Gaussian fits. 

As an additional test of the interpretation of the Gaussian fit parameters, we can compare 
their values to limits placed on cold Hi fractions based on double Gaussian fitting to individual 



line-of-sight profiles. Warren et al. (2012) estimated the spatially-resolved minimum and maximum 
fraction of cold Hi, J-coid, characterized by individual line of sight profiles that are best fit by a 
double Gaussian whose narrow component width is cr < 6 km s~^, in a number of our sample 



galaxies. In Figure 34 we plot the limits on An/Af, from Warren et al. (2012) versus A„/Afe values 
measured from double Gaussian fits to the global superprofile. Each box represents a single galaxy. 
The position and width of the box on the x-axis is determined by the A„/Af, value and associated 
error from double Gaussian fits. The size of the box on the y-axis is determined by the allowed 



range of An/Ah = Tcoid/{^ — ^coid) values for that galaxy as given by Warren et al. (2012). Grey 



boxes represent measurements from galaxies presented in this paper, while red boxes are those from 



12012 For simplicity, we have assumed that our A„/Aft uncertainty is 0.05. 
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The dashed hne in Figure 34 is a hne of equahty, where the independent hmits on A„/Afe 
match double Gaussian A„/Af, values. That is, the boxes would roughly lie along the dashed line 
if double Gaussian A„/A;, values matched the independent limits. However, this behavior is not 
seen. Instead, the independent limits on A^ / A5 tend to be smaller than those measured by double 
Gaussian fits. While we might expect the independent limits on A„/Af, to be smaller than the 
double Gaussian A„/Ab values because the superprofiles are a higher S/N representation of Hi 
spectra, we would at least expect to see the a positive correlation between independent limits and 
double Gaussian An/Ab values. However, the data, especially the red boxes, show, if anything, a 
negative correlation. The fact that double Gaussian A„/Aft values do not match independent limits 
indicates that the cold and warm H i interpretation may not be valid. 



Finally, we can examine limits on the H i velocity dispersions from observational studies. Petric 



h Rupen (2007) found that double Gaussian fits to median Hi profile shapes at different radii 
exhibited the same ratio of narrow to broad components. This finding complicates the CNM/WNM 
interpretation, as Hi in the outskirts of disks would not be expected to have the same fraction of 
gas in the warm component as regions inside r25. Second, the CNM/WNM interpretation of A„/Afc 
values implicitly assumes that there are two well-thermalized H i components in the ISM. If this were 
the case, we might expect to see two characteristic velocity dispersions in the statistical ensemble 



of individual Hi spectra widths. Braun et al. (2009) measured the non-thermal velocity dispersion 



of Hi in M31 and found no evidence of a bimodal distribution. Instead, they find a range of 
non-thermal velocity dispersions between 3 — 25 km s~^, with most line-of-sight profiles exhibiting 
widths ~ 8 km s~^. Based on this evidence, it therefore is more physically meaningful to interpret 
the central peak of the superprofiles, which show similar velocity dispersions of ~ 8 km s~^, as 
turbulence in the ISM instead of as the composite of emission from distinct cold and warm Hi 
components. However, this last constraint may not apply if the CNM and WNM components of 
the ISM are well-mixed on scales smaller than the spatial resolution of the observations. 

With the strength of the evidence presented above, we have chosen to parameterize the super- 
profiles as composed primarily of a central turbulent peak with wings to either side, instead of as 
two Gaussian components. The parameters of double Gaussians, when interpreted as representa- 
tive of the CNM and WNM, do not behave as expected when exploring the superprofiles in more 
detail, as the fit is dominated by the small amount of H i in the wings. Additionally, evidence from 
other studies does not support the cold and warm gas interpretation of the double Gaussian fits. 
Previously-derived limits on A^/A;, do not match A„/Afe values from the CNM/WNM interpreta- 
tion of double Gaussian fits, and there is no evidence for a bimodal temperature distribution in 
individual H i line of sight spectra even though such an interpretation is relevant in localized re- 
gions. We therefore believe the superprofiles are better interpreted not as two Gaussian components 
representing the CNM and WNM, but instead as a central turbulent peak with wings generated by 
kinematically disturbed gas. 
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Fig. 1. — Correlations between the globally- averaged physical properties for our sample. Each panel 
shows the correlation between two different properties. Within each panel, each point represents 
the globally-averaged value for a single galaxy. Points are colored black if the two properties are 
significantly correlated {ps < 0.01, or rs > 0.53) or grey if they are uncorrelated. 
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Fig. 2. — Two observed line-of-sight spectra from standard cubes with Upeak methods: the first 
moment map, a Gaussian fit, and a Gauss-Hermite fit. The upper panels show the spectra and 
various fits, while the lower panels are the residuals. The thin black line in the spectrum, the 
dashed line is the Gaussian fit, and the thick black line is the Gauss-Hermite fit. The plotted Hi 
profiles are for line-of-sight spectra with higher than average S/N ratios, to better show the adopted 
functional form. The Gauss-Hermite polynomials are primarily to find Upeak and are not meant to 
characterize the detailed line profile structure. Left: Sextans A, for a line-of-sight spectrum with 
S/N = 16.9. Right: Holmberg I, for a line-of-sight spectrum with S/N = 9.7. 
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Fig. 3. — The fpeak map generated from the Gauss-Hermite fits with the systemic velocity removed 
(top) and differences from both a single Gaussian fit (middle) and the first moment map (bottom). 
The red contour line shows regions with S/N> 5. 
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Fig. 4. — Example Gauss- Hermite polynomial fits to various Hi line-of-sight spectra for NGC 2366, 
a galaxy with a velocity resolution of 2.6 km s~^. The line-of-sight profiles are sorted into rows 
based on their S/N; spectra with S/N< 5 have grey backgrounds. The general asymmetry of the 
line-of-sight spectra is readily apparent. For spectra with S/N> 5, the Gauss-Hermite polynomials 
generally do a good job at finding the peak. At lower S/N, the peak is more difficult to determine 
and may even be due solely to noise spikes. In our analysis, we only use pixels with S/N> 5. 
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Fig. 5. — Same as Figure [4| for Sextans A, a galaxy with a velocity resolution of 1.3 km s ^. 
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Fig. 6. — The average error on finding t^pcak (o'tipcak) ^ function of S/N for three different simulated 
galaxy data cubes with three different velocity resolutions. The grey box indicates regions where 
'^'^peak is below the observation's velocity resolution. The dashed black line is our adopted S/N = 
5 threshold. At our S/N threshold, we find an average uncertainty of < 1 km s"^ for Av = 0.6 
km s-i (GR 8), < 1 for Av = 1.3 km s-^ (Sextans A), and < 1.5 for Av = 2.6 km s-^ (UGC 4483). 
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Fig. 7. — Hi superprofiles. In each panel, the black line represents the measured superprofile 
for each galaxy. Grey regions around the superprofile are the 1-a uncertainties on the flux. The 
dashed red line shows a Gaussian scaled to the amplitude and the half-width at half-maximum of 
the superprofile. Shaded red regions between this line and the superprofile are the "wing" regions 
and represent /wings- The vertical red line is the characteristic velocity of the wings, Swings • As the 
superprofiles are the analogue of integrated Hi spectra, but with the rotational velocity removed, 
we plot flux in Jy versus offset velocity. However, the Jy value is not indicative of our signal in a 
single channel. Galaxies are ordered by decreasing Mbaryon.tot- 
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Fig. 7. — Hi superprofiles (continued). Galaxies are ordered by decreasing Mb, 
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Fig. 7. — Hi superprofiles (continued). Galaxies are ordered by decreasing Mb, 



Offset Velocity / HWHM 



Fig. 8. — Each solid black line is the superprofile for a single galaxy. They have been normalized to 
their maximum height on the y-axis. The velocity axis has been scaled by the superprofile HWHM 
value. The superprofile lines have been plotted with transparency; overlap regions are darker. The 
thick dashed line is a Gaussian with amplitude = 1 and HWHM = 1. The lower panel shows the 
residuals from the Gaussian overlay. Compared to the Gaussian, the superprofiles are slightly more 
narrowly-peaked and have more flux in the wings. The superprofiles show a remarkably uniform 
shape, with the primary variations in the amount of gas moving faster than expected from the 
Gaussian overlay. 




Fig. 9. — The distribution of measured parameters (Tcentrab "Swings, /wingsi and a in our sample 
galaxies. 
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Fig. 10. — A comparison of the global second moment versus the dispersion measured by the 
superprofile HWHM-scaled Gaussian. The dashed line shows where the = (Tcentrai- The error 
bars on the x-axis are the approximate uncertainties on i7centraij while those on y-axis represent the 
weighted standard deviation of second moment values in the pixels considered. The solid black line 
at CTina /^central ~ 1-2 On the lower panel indicates the median ratio between Um2 and (Tcentrai- The 
second moment overestimates the width by ~ 10 — 50%. 
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Fig. 11. — The upper panel shows the measured parameters versus galaxy Mbaryon.tot- Each point 
is colored from blue to red based on increasing Mbaryon.tot- The bottom panel shows the residuals 
of the normalized superprofiles. The color of the line corresponds to that on the upper panel. The 
lines are plotted with a transparency value, so overlap regions are more saturated. Both /wings and 
o show trends with increasing star formation. 
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Fig. 12. — Same as Figure 



11 



but for SFR. We find that a decreases with increasing SFR. 



- 72 - 




Fig. 13. — Same as Figure 11, but for (Ssfr)- We find that galaxies with higher (Ssfr) have 
higher characteristic wing velocities. 
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Fig. 14. — Same as Figure 11, but for SFR / Mhi . Galaxies with higher SFR / Mhi also have a 
higher fraction of gas in the wings of their superprofiles. 
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Fig. 15. — The superprofile parameters as a function of (Shi)- We find that both o"centrai and Cwings 
increase with increasing (Shi)- 
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Fig. 16. — The superprofile parameters as a function of galaxy inclination. The measured parame- 
ters do not change systematically as a function of inclination, but we note that the inclinations for 
many of the dwarfs in our sample are very uncertain. 
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Fig. 17. — Observed superprofile properties versus physical properties when higher mass galaxies are 
included. The main dwarf sample is shown with filled black circles, while the higher mass galaxies 
are shown with open black circles. In most cases, the higher mass galaxies fit the correlations, but 
in a few, they do not (i.e., cTcentrai versus (Shi) and cJwings versus (Shi)). 
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Fig. 18. — Example galaxies in our sample that show asymmetric Hi line-of-sight profiles near 
star forming regions. The background image is Ssfr, with the black line representing the S/N > 5 
threshold. The color overlays represent asymmetric line-of-sight profiles, with blue indicating line- 
of-sight spectra where the first moment is smaller than Vpeak and red showing where the first 
moment is larger than Vpeak- The color ranges from -5 to 5 km s~^ differences. Regions with small 
absolute differences between Vpeak ^ind the first moment are also shown with more transparency, 
and those with large absolute differences are less transparent. We note that not all star forming 
regions are associated with asymmetrical line-of-sight profiles, and not all asymmetric line-of-sight 
profiles are near star forming regions. However, the observed overlap between some star forming 
regions and the strongest H l line-of-sight asymmetries may indicate that star formation can be one 
driver of asymmetry in Hi line profiles. 
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Fig. 19. — The left panel shows the energy available from the gravitational instability in Wada 



et al. ( 2002 ) over one turbulent timescale versus the turbulent energy in the central H i component 
as determined by the HWHM Gaussian fit (left). The entire background has been shaded in grey 
to represent the fact that all required efficiencies are unphysical, i.e., > 1. The right panel shows 
the associated implied efficiencies (e = -E^hi, central /-E'grav) versus input gravitational energy. It is 
clear that this instability cannot provide enough energy to drive the central peak in these galaxies 
over a single turbulent dissipation timescale. 
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Fig. 20. — The left panel shows the energy available from star formation over one turbulent 

timescale versus the central Hi component as determined by the HWHM Gaussian fit (left). The 
dashed grey lines indicate constant efficiency, and the shaded gray regions represent esF,turb > 0.1 
(light grey) and esF,turb > 1 (dark grey). The right panel shows esF,turb versus E^f- The mean 
efficiency required is esF,turb = 0.117 ± 0.096. 
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Fig. 21. — The left panel shows the energy available from star formation over 32.5 Myr versus Hi 
energy in the wings of the superprofile. The dashed grey lines indicate constant efficiency, and the 
shaded gray regions represent esF,wings > 0.1 (light grey) and esF,wings > 1 (dark grey). The right 
panel shows esF,wings versus £^SF- The mean efficiency required is esF,wings = 0.042 it 0.020. 
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Fig. 22. — The left panel shows the energy available from star formation over one turbulent 
timescale versus H I energy in the wings of the superprofile. The dashed grey lines indicate constant 

efficiency, and the shaded gray regions represent esF,wings > 0.1 (light grey) and esF,wings > 1 (dark 
grey). The right panel shows esF,wings versus E^f- The mean efficiency required is esF,wings = 
0.174 ±0.095. 
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Fig. 23. — The left panel shows the energy available from star formation over a single turbulent 
timescale versus Hi energy in the full superprofile. The dashed grey lines indicate constant effi- 
ciency, and the shaded gray regions represent esF,fuii > 0.1 (light grey) and esF,fuii > 1 (dark grey). 
The right panel shows esF.fuU versus £^SF- The mean efficiency required is esF,fuii = 0.275 it 0.190. 
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Fig. 24. — The left panel shows the energy available from star formation over 32.5 Myr versus Hi 
energy in the full supcrprofile component as determined by the HWHM Gaussian fit. The dashed 
grey lines indicate constant efficiency, and the shaded gray regions represent esF,full > 0.1 (light 
grey) and egF.fuii > 1 (dark grey). The right panel shows esF,fuii versus Esf- The mean efficiency 
required is esF,fuii = 0.102 ± 0.058. 
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Fig. 25. — A comparison of noise estimates for the superprofiles generated from the standard 
cube (black) and from the rescaled cube (red) for the galaxy GR 8. The solid lines represent the 
superprofiles, while the shaded regions show the 1-a noise estimate. The upper panels show the 
absolute flux measured in each superprofile, while the lower panels show the superprofiles after 
normalizing to the same maximum amplitude. The left panels shows the noise estimate given by 
Equation Al for both the standard and rescaled superprofiles. The right panels shows the noise 
estimate given by Equation Al for the standard superprofile and that given by Equation |8] for 
the rescaled superprofile, where the noise is scaled by the ratio of fluxes between the rescaled and 
standard superprofile. The rescaled noise approximately matches the fractional uncertainty of noise 
in the standard superprofile. 
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Fig. 26. — Fake superprofiles with included uncertainties in our Vpcak measurements. Each panel 
shows a different galaxy. The black line represents the fake superprofile. The dashed blue line is 
the input Gaussian and the dashed red line is the HWHM-Gaussian fit to the fake superprofile. 
All lines are plotted in the top panel, but the differences are smaller than the line widths. The 
bottom panel shows the residuals (i.e., observed - fit and observed - input). Uncertainties in Vpcak 
broaden the profile slightly, but the effects are ~ 100 times smaller than the observed amplitude of 
the wings. The effect due to uncertainties in fpeak is therefore negligible. Measurement errors in 
Upeak therefore do not create the shape of the observed superprofiles. 
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Fig. 27. — The resulting superprofiles for three NGC 2366 models without spatial smoothing 
applied (thick black line) and smoothed to 200 pc (thick dashed line). The observed superprofile 
for NGC 2366 is shown in grey. The upper panel is for the model using the observed rotation curve 
parameters and fixed velocity dispersion. The middle panel shows a rotation curve with a more 
extreme rise in the center, but is the larger amount of beam smearing produces no noticeable effect 
on the profile. The bottom panel shows a model with velocity dispersion that declines with radius 
based on the second moment map. In all cases, the effects of beam smearing are not strong enough 
to produce the wings or to substantially widen the intrinsic superprofile. 
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Fig. 28. — Superprofiles generated from binned cubes for DDO 125. The thick black hne represents 
the original superprofile. Colored lines represent superprofiles generated from cubes binned to 
Av = 1.3,2.6, and 5.2 km s~^. While the differences from the original superprofile are small for 
the Av = 1.3 and 2.6 km s~^ observations, the superprofile generated from the Av = 5.2 km s~^ 
cube is often different by more than 5% compared to the original. In the Av = 5.2 km s~^ cubes, 
the superprofiles are noticeably wider and shorter, thus leading to artificially inflated o"ccntrai values 
and decreased /wings values. 
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Fig. 29. — The change in each parameter (cxcentrah Cwingsj /wings, and a) with increasing velocity 
resolution. Each colored line represents a single galaxy. The parameters generally show only modest 
changes with velocity resolution. The most variable is a, whose variation has been accounted for 
in the uncertainties. 
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Fig. 30. — Scaled superprofiles for all the galaxies. The top panels scale the velocity axis by 
^central) the middle panels by (T„, and the bottom panel by ai,. The scaled velocity axes are then 
normalized such that the median scaled superprofile has the same width in all panels to better 
show the variation in shape. The median superprofile is shown in red. On the left we show the 
normalized superprofiles, and on the right we show the differences from the median superprofile. 
The HWHM scaling provides the best overall description of the shape of the superprofiles. 
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Offset Velocity / HWHM Offset Velocity / HWHM 

Fig. 31. — Double Gaussian fit shapes compared to the overall superprofile shapes. The black lines 
are the full superprofile, while the red and blue lines represent that narrow and broad Gaussian 
components of a double Gaussian fit. Both the superprofiles and Gaussian components have been 
scaled to the HWHM of that specific superprofile. The left panel shows the full sample, while 
the right panel highlights NGC 7793 and Sextans B to show the drastic difference in Gaussian 
components for superprofiles whose wing fluxes differ by only 9% when scaled to the same HWHM. 
This small change drives large differences in the best-fit double Gaussian components. 
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Fig. 32. — The relationship between Narrow / Broad Amplitude of double Gaussian fits compared 
to the width of the wings relative to the characteristic superprofile width measured by the HWHM. 
The grey points in the left panel are galaxies from our sample, while those in right panel (black) 
are taken from 12012 In both cases, there is a clear trend that galaxies with broader wings have 
relatively lower broad amplitudes. 
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Fig. 33. — The relationship between An/Af, compared to the fraction of gas in the wings, /wings- 
For galaxies in the 12012 sample, we find that galaxies with more gas in the wings have higher 
fraction of gas in the narrow component. If the A„/Aft parameter were tracing the ratio between 
the mass of CNM to WNM, we would on average expect galaxies with more flux in the wings of the 
profile to have a more Hi in the WNM, and therefore smaller A^/Afe ratios. This expectation is the 
opposite of how the A„/Aft parameter behaves. It is instead likely that A„/Ab, like the amplitudes 
of the double Gaussian components, is driven by the relative broadness of the wings. 
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Fig. 34. — A comparison between A^/Aft values from double Gaussian fits and independent limits 



on An/Aj, from Warren et al. (2012). Each box represents an individual galaxy. The x-position 
and width are determined by double Gaussian A„/Ab values and associated errors. The top and 



bottom of each box is the limit placed on A„/Afe by Warren et al. (2012). The line of equality is 
shown as a thick dashed line. The A^/A^ values measured by double Gaussian fits do not match 



the limits placed by Warren et al. (2012). 
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Table 2. Observation Parameters 



Of! law 


Av 


'^2UUpc 






(km s 


\ ) 


I'tti Tv Hpam ^\ 


1 


2 


3 


4 


NGC 7793 


2.6 


10.58 


1.19 


IC 2574 


2.6 


10.88 


0.91 


NGC 4214 


1.3 


13.57 


1.08 


Ho II 


2.6 


12.21 


1.54 


NGC 2366 


2.6 


12.85 


1.00 


DDO 154 


2.6 


9.59 


0.60 


Ho I 


2.6 


10.58 


1.23 


NGC 4190 


1.3 


11.79 


1.50 


NGC 3741 


1.3 


12.73 


1.90 


Sextans A 


1.3 


29.89 


4.06 


DDO 53 


2.6 


11.43 


0.76 


DDO 190 


2.6 


14.69 


0.62 


DDO 125 


0.6 


15.99 


3.07 


Sextans B 


1.3 


29.68 


1.90 


DDO 99 


1.3 


15.93 


1.97 


M81 DwB 


2.6 


7.78 


0.64 


UGCA 292 


0.6 


11.40 


2.13 


NGC 4163 


0.6 


14.42 


2.16 


UGC 4483 


2.6 


12.10 


0.80 


DDO 181 


1.3 


13.14 


1.52 


UGC 8833 


2.6 


13.39 


0.59 


DDO 187 


1.3 


18.67 


2.12 


GR8 


0.6 


19.83 


4.15 



Note. — (1) Galaxy name. (2) Channel spacing. 

(3) Circular beam in " that corresponds to 200 pc. 

(4) rms noise in mJy beam"-*^ in convolved, standard 
cube. 
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Table 4. Measured Superprofile Parameters 



CT?^l^^">^■v 


\^m2/ 
^km s 


'^central 

(^km S 


Swings 


Jw 




Nk 
beams 


1 


2 


3 


4 


5 


6 


7 


NGC 7793 


11.3 


8.1 ±0.5 


26.0 ± 1.4 


0.15 ±0.02 


0.04 ±0.06 


909 


IC 2574 


9.3 


7.2 ±0.5 


19.8 ± 1.0 


0.12 ±0.02 


0.11 ±0.06 


2443 


NGC 4214 


9.3 


6.5 ±0.5 


21.8 ± 1.2 


0.15 ±0.02 


0.05 ± 0.06 


1558 


Ho II 


9.0 


7.1 ±0.5 


21.3 ± 1.0 


0.12 ±0.02 


0.14 ±0.06 


925 


NGC 2366 


13.2 


10.1 ±0.5 


30.2 ±0.9 


0.13 ±0.02 


0.22 ± 0.06 


892 


DDO 154 


8.6 


7.5 ±0.5 


19.8 ±0.4 


0.09 ± 0.02 


0.09 ± 0.06 


682 


Ho I 


9.0 


6.7 ±0.5 


20.6 ±2.1 


0.15 ±0.03 


0.11 ±0.06 


162 


NGC 4190 


11.6 


9.3 ±0.6 


28.0 ± 7.2 


0.12 ±0.03 


0.12 ±0.07 


45 


NGC 3741 


8.5 


7.2 ±0.5 


20.3 ±3.0 


0.12 ±0.03 


0.15 ±0.07 


112 


Sextans A 


9.5 


8.4 ±0.5 


23.4 ± 1.0 


0.10 ±0.02 


0.09 ± 0.06 


215 


DDO 53 


10.3 


8.2 ±0.5 


23.8 ± 1.6 


0.11 ±0.03 


0.18 ±0.07 


93 


DDO 190 


10.4 


8.9 ±0.5 


24.6 ±3.1 


0.09 ±0.02 


0.25 ±0.08 


76 


DDO 125 


7.0 


6.0 ±0.5 


16.9 ±2.3 


0.14 ±0.03 


0.10 ±0.07 


53 


Sextans B 


7.9 


7.3 ±0.5 


17.8 ±0.7 


0.08 ±0.02 


0.26 ±0.07 


259 


DDO 99 


8.0 


7.3 ±0.5 


18.9 ± 1.6 


0.10 ±0.03 


0.19 ±0.08 


92 


M81 DwB 


12.5 


9.2 ±0.6 


30.7 ±5.0 


0.11 ±0.05 


0.23 ± 0.08 


20 


UGCA 292 


8.4 


7.8 ±0.5 


20.3 ±4.8 


0.06 ± 0.03 


0.27 ±0.10 


45 


NGC 4163 


8.5 


7.7 ±0.6 


22.1 ±3.1 


0.09 ±0.04 


0.40 ±0.12 


14 


UGC 4483 


9.9 


8.5 ±0.6 


24.9 ±2.6 


0.10 ±0.03 


0.41 ±0.11 


55 


DDO 181 


7.8 


6.5 ±0.5 


19.4 ±2.8 


0.14 ±0.04 


0.20 ±0.08 


54 


UGC 8833 


9.6 


8.0 ±0.6 


23.2 ±3.2 


0.11 ±0.03 


0.46 ±0.13 


30 


DDO 187 


11.3 


10.4 ±0.6 


25.0 ±3.3 


0.10 ±0.03 


0.22 ± 0.08 


23 


GR8 


8.0 


7.5 ±0.5 


20.2 ±2.4 


0.12 ±0.04 


0.29 ±0.10 


17 



Note. — (1) Galaxy name. (2) Average intensity- weighted global second moment 
value. (3) Width of central superprofile peak. (4) Characteristic velocity of the wings. 
(5) Fraction of H I in the wings. (6) Wing asymmetry parameter. (7) Number of inde- 
pendent beams contributing to superprofile. 
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Table 5. Double Gaussian Fit Parameters 



Name (7„ ai, An /At, Onjo^ 

(km s"-*^) (km s~^) 
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u. 
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Nnn 41 Qo 


7 
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-1- 


9 9 
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4- 

ZL 





1 7 





48 


-1- 

ZL 





04 


NGC 3741 


4. 


,7 


± 


0, 


.4 


10.5 


± 


0.8 


0, 


.34 


± 


0, 


.09 


0, 


.45 


± 


0, 


.01 


Sextans A 


6. 


,2 




0, 


,2 


12.3 




0.4 





.54 




0, 


.07 





.51 




0, 


.01 


DDO 53 


6. 


,8 




0, 


.1 


13.9 




1.1 


1, 


.01 




0, 


.10 





.49 




0, 


.03 


DDO 190 


6. 


,8 




0, 


.4 
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± 
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.53 




0, 


.19 





.55 
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.01 
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4. 


,2 




0, 
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.06 
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13.3 


± 


3.7 


1 


.32 


± 


0, 


.01 





.50 


± 


0, 


.10 


UGC 4483 


7. 


,0 




0, 


.2 


14.0 




1.4 


0, 


.98 




0, 


.12 


0, 


.50 




0, 


.03 


DDO 181 


5. 


,0 




0, 


.7 


10.3 




2.0 


0, 


.60 




0, 


.42 


0, 


.48 




0, 


.01 


UGC 8833 


6. 


,7 




0, 


.2 


13.2 




1.5 


1 


.01 




0, 


.11 





.50 




0, 


.05 


DDO 187 


5. 


,3 




0, 


.3 


13.4 




0.7 





.15 




0, 


.03 





.40 




0, 


.01 


GR8 


5. 


.3 


± 


0, 


.6 


11.1 


± 


1.6 


0, 


.44 


± 


0, 


.19 


0, 


.48 


± 


0, 


.02 



Note. — All double Gaussian parameter uncertainties calcu- 
lated from noise uncertainties only.(l) Galaxy name. (2) Width 

of narrow Gaussian component. (3) Width of wide Gaussian com- 
ponent. (4) Ratio of area of narrow Gaussian component to broad 
Gaussian component. (5) Ratio of narrow width to broad width. 
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Table 6. Spearman Correlation Coefficients 



<^central Cwings Jwings 



Property 


rs 


PS 


rs 


PS 


rs 


PS 


rs 


PS 




0.079 


0.720 


0.261 


0.229 


0.262 


0.227 


-0.459 


0.027 


l^baryonjtot 


-0.189 


0.388 


0.083 


0.707 


0.491 


0.017 


-0.769 


<0.001 




-0.150 


0.494 


0.085 


0.700 


0.414 


0.050 


-0.655 


<0.001 


Mhi 


-0.152 


0.488 


0.047 


0.830 


0.377 


0.076 


-0.643 


<0.001 


SFR 


-0.071 


0.747 


0.221 


0.310 


0.486 


0.019 


-0.721 


<0.001 


SFR / Mhi 


0.066 


0.764 


0.390 


0.066 


0.529 


0.010 


-0.179 


0.414 




0.297 


0.168 


0.626 


0.001 


0.314 


0.144 


0.000 


1.000 


(Shi) 


0.536 


0.008 


0.707 


<0.001 


-0.053 


0.809 


0.200 


0.361 


(^baryon) 


0.336 


0.117 


0.623 


0.002 


0.123 


0.578 


0.115 


0.603 


i 


-0.218 


0.318 


-0.311 


0.149 


0.122 


0.579 


-0.278 


0.199 



Note. — Spearman correlation coefficient and probability ps between superprofile 
parameters and physical properties. Significant correlations (i.e., ps < 0.01) are shown 
in bold. 
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Table 7. Measured Superprofile Parameters for More Massive Galaxies 



Galaxy 




'-'^central 


'Swings 


fw 


a 


■^beams 




(km s^"*^) 


(km S^^) 










1 


2 


3 


4 


5 


6 


7 


NGC 5055 


13.6 


8.8 


34.3 


0.18 


0.05 


3616 


NGC 2903 


16.6 


8.8 


39.5 


0.21 


0.06 


2186 


NGC 5236 


10.1 


8.8 


29.0 


0.12 


0.18 


447 


NGC 3351 


13.5 


6.6 


29.0 


0.17 


0.05 


607 


NGC 4736 


14.8 


8.6 


40.6 


0.20 


0.04 


725 


NGC 628 


7.5 


6.5 


18.3 


0.08 


0.12 


2819 


NGC 2403 


10.0 


8.3 


23.0 


0.09 


0.05 


2253 


NGC 2976 


14.5 


12.3 


36.0 


0.10 


0.18 


110 



Note. — Measured superprofile parameters for higher-mass 
THINGS spirals. Columns are the same as Table [4j 
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Table 8. Implied Scale Heights 



VTCtid-A-Y 


\'I'Z/ 




(no) 


1 

± 





NGC 7793 


78 ±43 


IC 2574 


336 ± 87 


NGC 4214 


228 ± 103 


Ho II 


302 ± 75 


NGC 2366 


479 ± 93 


DDO 154 


708 ± 139 


Ho I 


191 ± 41 


NGC 4190 


130 ± 65 


NGC 3741 


561 ± 123 


Sextans A 


438 ± 84 


DDO 53 


259 ± 51 


DDO 190 


214 ±65 


DDO 125 


329 ± 135 


Sextans B 


639 ± 224 


DDO 99 


617 ± 161 


M81 DwB 


140 ± 65 


UGCA 292 


172 ± 29 


NGC 4163 


120 ± 61 


UGC 4483 


322 ± 67 


DDO 181 


346 ± 104 


UGC 8833 


313 ± 76 


DDO 187 


514 ± 108 


GR 8 


344 ± 100 



Note. — Implied scale 
heights for sample galax- 
ies. (1) Galaxy name. 
(2) Average scale height 
implied by (Tcentrai and 
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(Shi). 



